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EXPLICIT VON NEUMANN STABILITY CONDITIONS FOR 
THE c-r SCHEME — A BASIC SCHEME IN THE DEVELOPMENT 
OF THE CE-SE COURANT NUMBER INSENSITIVE SCHEMES 


Sin-Chung Chang 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Abstract 


As part of the continuous development of the space-time conservation element and solution element 
(CE-SE) method, recently a set of so called “Courant number insensitive schemes” has been proposed. The 
key advantage of these new schemes is that the numerical dissipation associated with them generally does 
not increase as the Courant number decreases. As such, they can be applied to problems with large Courant 
number disparities (such as what commonly occurs in Navier-Stokes problems) without incurring excessive 
numerical dissipation. 

A basic scheme in the development of the Courant number insensitive schemes is the so called “c-r 
scheme” . It is a solver of the PDE 


du du 

m +a d^ = 0 

where 0 is a constant. At each space-time staggered mesh points (j, n), the c-r scheme is formed by 


( U ^d/2 ~ ( U ^+l/l]} 


and 


(«*)* 


1 

2(1+7) 


n— 1/2 
“j+1/2 


(1 + 2 " - T ){ u x) n j+ l/l 


— u 


n- 1/2 

j i/2 


( 1 - 2 V- t){Ux) 


n— 1/2 
3- 1/2. 


Here: (i) u ™ and (i^)™, respectively, denote the numerical analogues of u and (ax/4 )du/dx at the mesh 

point (j, n); (ii) v = f aAt / ax is the Courant number; and (iii) r is an adjustable parameter 7 ^ —1. 

Because the c-r scheme is formed by two rather complicated equations involving two parameters v and 
r, it were not expected that its von Neumann stability conditions could be cast into an explicit analytical 
form. Against this expectation, it will be shown rigorously in this paper that, based on the von Neumann 
analysis, the c-r scheme is stable if and only if 

V 1 < 1, t>To(v 2 ), and (zv 2 ,t)^(1,1) 


where 


ro 


MX) d = f < 


if x = 0 


4 — x — 2 ^ 2(2 — x — x 2 ) 


- 1 + Vl - 2x + 5x 2 
2x 


if 0 < x < 3/11 


if 3/11 <x<l 


Note that the current stability conditions are in complete agreement with those generated numerically and 
reported earlier. 

In addition, it will be shown that: (i) r G (x) is continuous at x = 0; (ii) r G (x) is consistently defined at 
x = 3/11; (iii) 

lim r'(x) = lim r'(x) = 121/90 


where r'(x) ^= f dr 0 {x)/dx ; (iv) r G (x) is strictly montonically increasing in the interval 0 < x < 1 ; and (v) 

x < r 0 (x) < \fx , 0 < x < 1 
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1. Introduction 


As part of the continuous development of the space-time conservation element and 
solution element (CE-SE) method [1-11], recently a set of so called “Courant number 
insensitive schemes” has been reported in [9-11]. The key advantage of these new schemes 
is that the numerical dissipation associated with them generally does not increase as the 
Courant number decreases. As such, they can be applied to problems with large Courant 
number disparities (such as what commonly occurs in Navier-Stokes problems) without 
incurring excessive numerical dissipation. 

A basic scheme in the development of the Courant number insensitive schemes is the 
so called “c-r scheme” [11]. It is a solver of the PDE 


du du 

at +a te =0 


where a ^ 0 is a constant. Consider Fig. 1 and let Q denote the set of all space-time 
staggered mesh points (dots in Fig. 1), where n = 0, ±1/2, ±1, ±3/2, ±2, . . ., and, for each 
n, j = n ± 1/2, n ± 3/2, n ± 5/2, . . .. Then, at each (j, n ) G 0, the c-r scheme is formed by 


u j 2 


{(1 + U ) ul j - 1/2 + (1 “ U ) U ' j + l /2 + (1 “ y2 ) (Ab)J_i/2 “ ( M *) j + l /2 } 


and 


2(\ + t) [“i+1/2 - (! + 2l/ - r )(“i)"+i 1 /2 2 - “"-1/2 - (! - 2 “ - 

(1.3) 


Here: (i) u™ and (ux)'/, respectively, denote the numerical analogues of u and (ax/ 4 )du/dx 
at the mesh point (j, n ) ; (ii) 


def CL At 
V = 

AX 


(1.4) 


is the Courant number; and (iii) r is an adjustable parameter ^ — 1. It is shown in [12] 
that Eqs. (1.2) and (1.3) are consistent with a pair of PDEs with Eq. (1.1) being one of 
them. 


Because the c-r scheme is formed by two rather complicated equations involving two 
parameters v and r, it was not expected that its von Neumann stability conditions could 
be cast into an explicit analytical form. But to the contrary, it will be shown rigorously in 
this paper that, based on the von Neumann analysis, the c-r scheme is stable if and only 
if 

v 2 < 1, t > t 0 (u 2 ), and (z/ 2 , r) ^ (1, 1) (1.5) 


where 



x = 0 

0 < x < 3/11 
3/11 < x < 1 


( 1 . 6 ) 
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Note that the current stability conditions are in complete agreement with those generated 
numerically and reported earlier in [11]. 

In addition, it will be shown that: (i) r 0 (x) is continuous at x = 0; (ii) t 0 (x) is 
consistently defined at x = 3/11; (iii) 

lim t' 0 (x) = lim t' 0 (x ) = 121/90 (1-7) 

X ^TT X ^TT 

where (a:) dr 0 (x)/dx; (iv) t 0 (x) is strictly monotonically increasing in the interval 
0 < x < 1; and (v) 

x < r 0 (x) < \/x, 0 < x < 1 (1.8) 

Eqs. (1.5) and (1.8) coupled with the facts that r o (0) = 0 and \fv ^ = p| imply that 
the c-r scheme is stable if 

T = \u\ < 1 (1-9) 

On the other hand, Eqs. (1.5) and (1.8) imply that the c-r scheme is unstable for the cases 

i) 

v 2 > 1 (1.10) 

and (ii) 

r = v 2 and 0 < v 2 < 1 (1.11) 

Note that, for a reason explained in [9,11], the special c-r scheme with Eq. (1.9) is a 
Courant number insensitive solver for Eq. (1.1). 

The rest of the paper is outlined as follows. For any pair of v and r, and any phase 
angle 9, the amplification matrix Q ( u. r, 9) that arises from the von Neumann stability 
analysis is presented in Sec. 2 (see Eq. (2.8)). The definition of stability (Definition 1) 
is then given in the same section in terms of the behaviors of [Q(is,T,9)] m , —7 r < 9 < n, 
as the integer m — * +oo. In Sec. 3, Theorems 1 and 2 are introduced to link stability 
with the spectal radii p(Q(u,r,9)) of Q(u,r,9), —it < 9 < it. Based on the preliminaries 
given in Secs. 2 and 3, the main results are given in Sec. 4. Specifically, Sec. 4 begins 

with Theorem 3, in which the necessary and sufficient stability conditions are expressed 

implicitly in terms of a requirement on p(Q(i/,r,9)), —tt < 9 < tt. It is then followed 
by a systematic and rigorous effort to obtain the explicit solution to the above implicit 
conditions. Finally, conclusions and discussions are presented in Sec. 5. Moreover, to give 
the reader extra confidence on the main results established analytically in Theorems 34 and 
35, these theorems are further validated numerically in Appendices A and B, respectively. 
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2. von Neumann Stability Analysis 


For any (j. n ) € fl, let 


and 


where 


->( • \ def 


u. 


Ml 


^ / n def 1 

Q+wv = - 


r) t n def 1 

Q-fav = - 


/l + U 

-1 


V 1 + r 

(1-v 
1 


Vi + 


1-u 2 \ 

1-2 V-T 

1 + T ) 

1+2 V — T 
1 + r / 


1 + r^O 

is assumed. Then Eqs. (1.2) and (1.3) can be expreseed as 

q(j,n) = Q+q(j ~ 1/2, n - 1/2 ) + Q-q(j + 1/2, n - 1/2) 


( 2 . 1 ) 

(2.2) 

(2.3) 

(2.4) 

(2.5) 


Hereafter Q + (u,t) and Q-(v,r) may be abbreviated as Q + and Q-, respectively. 

To study the stability of the c-r scheme using the von Neumann analysis [1], for all 
(j, n ) G f2, let 

q(j,n) = q*(n,9)e lj6 (2.6) 

Here (i) i \[—l, (ii) 9, — oo < 9 < +oo, is the phase angle variation per ax, and (iii) 
q*(n , 9) is a 2 x 1 column matrix. Substituting Eq. (2.6) into Eq. (2.5) and using Eq. (2.4), 
one has 

q*(n + l/2,6) = Q(v,T,6)q*(n,9) (2.7) 

where n = 0, ±1/2, ±1, ±3/2, . . ., and 


Q(v, t,9) = f e~ t6/2 Q + (p, t) + e ld/2 Q-(u,r) 

/ cos(#/2) — iusm(9 /2) 

~~ isin(0/2) 

V 1 + T 

Because of Eq. (2.7), Q{y , r, 9) is referred to as the amplification matrix of the c-r scheme 
per marching step (or per At/2). Also, by using Eq. (2.7), one has 

q*(n + m/2, 9) = [Q(v, r, 9)] m q*(n, 9) (2.9) 


— i{l — v 2 ) sin(0/2) \ 

(1 — r) cos(#/2) + 2iusin(9/2) 


1 + r 


(2.8) 


) 
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where m = 1, 2, 3, . . . and n = 0, ±1/2, ±1, ±3/2, . . .. 

As a result of Eq. (2.9), we have Definition 1. 

Definition 1. The c-t scheme is said to be stable with respect to a given ordered pair 
(v, r) if, for every 9 , — oo < 9 < +oo, all elements of the matrix [Q(n, r, 9)] m associated 
with this pair remain bounded as the positive integer m — * +oo. On the other hand, the 
scheme is said to be unstable with respect to a given (u, r) if, for any 9. — oo < 9 < +oo, at 
least one element of the matrix [Q ( u, r, 9)] in associated with this {+. r) becomes unbounded 
as to — * +oo. Hereafter, a given (u, r) is said to be c-t stable (unstable) if the c-t scheme 
is stable (unstable) with respect to this ( v , t). 

Note that: (i) Eq. (2.8) implies that, for any integer £, 

Q(v , t, 9 + 2£ir) = (-1 YQ(v, t, 9) (2.10) 

and (ii) for any 9, — oo < 9 < ±oo, there are a 9 ' , — it < 9' < ir and an integer £ such 
that 9 = 9' + 2£n. As such, Definitions 1 is equivalent to the simplified form in which the 
original range of 9 , i.e., — oo < 9 < +oo, is replaced by 

—7 T < 9 < 7T (2-11) 

Hereafter, the simplified form of Definition 1 is assumed. 

Given Definition 1, it will be shown in this paper that a given (u, r) is c-t stable if and 
only if it satisfies Eq. (1.5). As a first step, in Sec. 3 we will answer the following question: 
For any given ordered set (u,t,9), what are the requirements the matrix Q(u,t,9) must 
meet so that all elements of the matrix [ Q(v , r, 9)] m will remain bounded as to — » +oo? 
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3. Two Matrix Theorems 


Let M be any N x N matrix with real or complex elements. By definition, the 
eigenspace of M is the vector space spanned by its eigenvectors. Let the dimension of this 
eigenspace be denoted by N' . Then 1 < N' < N . The matrix is said to be (i) nondefective 
if N' = N and (ii) defective if N' < N [13]. 

Hereafter let N = 2. Then the eigenvalues Ai and A2 of the matrix M are the two 
roots of a quadratic characteristic equation. Moreover, we have Theorem 1. 

Theorem 1. The matrix M is defective if and only if (i) Ai = A2, and (ii) M 7^ A C I, 
where I is the 2x2 identity matrix and A c is the common value of Ai and A2. 

Proof. Let b\ and 62 be two nonnull 2x1 column matrices with 

Mb t = X e b e , £ = 1,2 (3.1) 

— > 

Then, for each &, b( is an eigenvector of M with the eigenvalue A/ : . In case that Ai 7^ A2, it 
— > — ^ 

is known that b 1 and 62 are linearly independent [13]. Thus N' = 2 and M is nondefective. 

Next let Ai = A2 and M be nondefective. Then N' = 2, i.e., there exist two linearly 

— ^ — > 

independent 2x1 column matrices b\ and 62 that satisfy Eq. (3.1). Let 



£ = 1,2 


and 



Then, because Ai = A2, Eq. (3.1) can be expressed as 


(Af - A C I)B = 0 


(3,2) 


(3,3) 


(3.4) 


where A c is the common value of Ai and A2. Because b\ and 62 are linearly independent, B 
is nonsingular [13]. Thus, B~ x , the inverse of B. must exist. Multiplying the expressions on 
the two sides of Eq. (3.4) from the right with B~ l leads to the conclusion that M — \ C I = 0, 
i.e., M = A C I . 

Conversely let M = A C I where A c is any scalar. Then it can be shown easily that 
(i) Ai = A2 = A c , and (ii) any 2x1 nonnull column matrix is an eigenvector of M. The 
conclusion (ii) implies that N' = 2 and thus M is nondefective. 

It has been shown that: (i) M is nondefective if Ai 7^ A2; and (ii) in case that Ai = A2, 
M is nondefective if and only if M = A C I (i.e., M is defective if and only if M ^ A C I) 
where A c is the common value of Ai and A2. Thus the proof is completed. QED. 
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Next let (i) m be an integer > 0; and (ii) p(M) be the spectral radius of M , i.e. , 


p(M) = max{|Ai|, |A 2 |} 


(3.5) 


Then we have Theorem 2. 

Theorem 2. Every element of M m will remain bounded asm-> +oo if and only if 


p(M) 


< 1 if M is nondefective 

< 1 if M is defective 


(3.6) 


Proof. According to the Jordan canonical form theorem [13], there exists a nonsin- 
gular 2x2 matrix S such that 

M = SAS- 1 (3.7) 

Here (i) 5 1-1 is the inverse of S ; (ii) 

/Ax 0 \ 

A = f I I if M is nondefective (3.8) 

\o a 2 y 

and (iii) 

A = f | J if M is defective (3.9) 

\0 A c ) 

Note that A c in Eq. (3.9) is the common value of Ax and A 2 in the defective case. 

By using Eqs. (3.8) and (3.9), one has: (i) 


/Af 0 \ 

A m =1 J if M is nondefective 

V 0 A ?) 


and (ii) 

(A™ mA™- 1 \ 

A m =| J if M is defective 

V 0 A™ / 

Because (i) Eq. (3.7) implies that 

M m = SA m S~ 1 


and (ii) Eq. (3.12) is equivalent to 


A m = 


(3.10) 


(3.11) 


(3.12) 


(3,13) 
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one can infer from Eq. (3.10) that, for the nondefective case, every element of M' n will 
remain bounded as m — » +oo if and only if 

p(M) < 1 (the nondefective case) (3-14) 

On the other hand, for the defective case, by using (i) p(M ) = |A C |, and (ii) 

TO if |A C | < 1 

lim |mA! n_1 | = < (3.15) 

m ^ + °° 1 +oo if |A C | > 1 

Eqs. (3.11)-(3.13) imply that, for the defective case, every element of M m will remain 
bounded as M — > +oo if and only if 

p(M) < 1 (the defective case) (3.16) 

Because Eq. (3.6) is the combined form of Eqs. (3.14) and (3.16), the proof is completed. 

QED. 

At this juncture, note that the term |mA™ _1 | grows linearly with rn as rn — +oo if 
| A c | = 1. Thus, for the defective case with |A C | = 1, the growth rate of the magnitude of 
any element of M m as m — > +oo is very low compared with the exponential growth rate 
associated with a nondefective or defective case with p(M) > 1. The implication of this 
observation will be addressed later. 
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4. Main Results 


An immediate result of Definition 1 and Theorem 2 is Theorem 3. 

Theorem 3. A given (u, r) is c-r stable if and only if the condition 

{ <1 if Q(u, r, 9) is nondefective 

(4T) 

<1 if Q(v, r, 9) is defective 

associated with the given (y, r) is met for all 9 , —7 r < 9 < n. 

Two immediate results of Theorem 3 are Theorems 4 and 5. 

Theorem 4. A necessary condition for any given (u. r) to be c-r stable is 

p{Q{y , r, 9)) < 1 , -7r < 9 < 7T (4.2) 


Theorem 5. In case that 

p(Q(v,r,9)) ± 1 (4.3) 

for all defective Q(u,t,9 ) (— r < 9 < n) associated with a given (zq r), Eq. (4.2) is also a 
sufficient condition for this (u, r) to be c-r stable. 

From Theorem 3, it becomes clear that a thorough stability study of the c-r scheme 
requires a systematic investigation of the matrix Q(v, r, 9) and its eigenvalues over the 
entire range of u, r, and 9. In the following, first we shall try to narrow down the possible 
(u,t) that are c-r stable by ruling out those that fail to satisfy Eq. (4.2). 

Let det(M) denote the determinant of any square matrix M . Then any eigenvalue A 
of Q(u, t, 9) satisfies the characteristic equation det (Q(u, r, 9) — XI) = 0, i.e., 

(1 + r)A 2 — [2r cos(#/2) — iu ( 3 + r) sin(0/2)] A 

— (1 — r) cos 2 (#/2) — (1 + i/ 2 ) sin 2 (6 , /2) — iv{l + r) sin(#/2) cos(#/2) = 0 
Let 

X(u, r, 9) = f 4 cos 2 (0/2) + [4(1 + r) — v 2 (t 2 + 2 r + 5)] sin 2 (0/2) (4.5) 

and 

Y (u, r, 9) ^= f 4^(1 — r) sin(0/2) cos(0/2) (4-6) 

Then, with the aid of Eq. (2.4), Eq. (4.4) implies that A = A + (u,t,9) or A = A _(v,t,9) 
where 


. , def 2tcos(9/2) - iu ( 3 + r) sin(0/2) ± \/X + W 

A±(i '’ T '" ) = 20TT) ! 


1 + T ^ 0 


(4.7) 
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Hereafter X(u, r, 9) and Y (//, r, 9) may be abbreviated as X and Y . respectively. Because 
the range of the phase angle 0 in the polar form of the principal square root \/X + iF is 
— 7 t/2 < 0 < 7t/ 2, it can be shown that 


VX + iF 


1 

75 


\[X ‘‘ 2 + F 2 + X + isign(F ) 7 \[Y 2 + F 2 - X 


where 


sign(F) 


I 

-1 


if F > 0 
if F < 0 


With the aid of Eq. (4.8), Eq. (4.7) implies that 


(4.8) 


(4.9) 


X±(u,t,9) = 


1 


2(1 + r) 


2 r cos(6>/2) ± V V^ 2 + F 2 + X 

v 2 


i i/(3 + r) sin(0/2) =F sign(F)^/ \/ X 2 + F 2 — X j> 


(1+r^O) (4.10) 


Next Eq (4.10) is used to yield 

2(1 + t) 2 (\A + j 2 + \A_I 2 ) = 4 t 2 cos 2 (9/2) + v 2 (3 + t) 2 sin 2 (9/2) + TX^+Y* (4.11) 


and 


(1 + r) 2 |A + | 2 |A_| 2 = (1 - r) 2 cos 4 (6»/2) + (1 + v 2 ) 2 sin 4 (6>/2) 


, 2\2 • 4 / 


(4.12) 


+ (2 — 2t + 3i/ 2 + t 2 v 2 ) sin 2 (#/2) cos 2 (#/2) 

For simplicity, hereafter A + (u,t,9) and A -(u,r,9) may be abbreviated as A + and A_, 
respectively. Next, let 

s *= f sin 2 (#/2), —7 r < 9 < n (4-13) 

Then 

cos 2 (6/2) = 1 — s (4-14) 

and, corresponding to the domain — ir < 9 < ir, the range of s is 


0 < s < 1 


(4.15) 


Next, let 


D(iy, r, s) = 2(1 — z/ 2 )(r 2 — v 2 )s 2 + [4 r + (r 2 — 6r — 3)z/ 2 ] s + 4, 0 < s < 1 

E(u, r, s) d = [16 r 2 - 8 (r 3 + 4r 2 + r + 2)v 2 + (r 2 + 2r + 5) V] s 2 
+ 8 [4t + (r 2 — 6r — 3 )u 2 ~\ s + 16, 0 < s < 1 


(4.16) 

(4.17) 
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and 


F(u, r, s) = (1 — v 2 ){y 2 — t 2 )s 2 — [2t(1 — t) + (3 + r 2 )z/ 2 ] s + 4 r, 0 < s < 1 (4.18) 

Then, by using Eqs. (4.5), (4.6), and (4.11)-(4.14), it can be shown that 

E(v,T,s) = [X(i ',T,0)f + [Y(v,t,0)] 2 > 0 (4.19) 

D(v, t, a) - ^/E(v,r,s) = 2(1 + rf (1 - |A + | 2 ) (1 - |A_ | 2 ) (4.20) 

and 

F(v,t, s) = (1 +r) 2 (1 - |A+| 2 |A-| 2 ) (4.21) 

As a preliminary to the future development, let 

H(y, r, s) d = [D(u, r, s)] 2 - E(u, r, s) (4.22) 

Then Eqs. (4.16) and (4.17) imply that 

H (v, r, s ) = 4(1 — p 2 )s 2 G(u, r, s) (4.23) 


where 

G(u, r, s) = f (1 - ^ 2 )(r 2 - v 2 ) 2 s 2 + (r 2 - v 2 ) [u 2 r 2 + (4 - 6 v 2 )t - 3 v 2 ] s 
+ 4 r \u 2 t 2 + (1 — v 2 )t — v 2 ] , 0 < s < 1 

With the above preparations, we have Theorem 6. 

Theorem 6. (A) For any (za r), the condition Eq. (4.2) is equvalent to the conditions 



D(v, r, s) > 0, 

0 < 5 < 1 

(4.25) 

H( is, r, s) > 0, 

0 < s < 1 

(4.26) 

F(v,t, s ) > 0, 

0 < s < 1 

(4.27) 


(B) Eqs. (4.25)-(4.27) are necessary conditions for any (u, r) to be c-r stable. 

Proof. Part B is an immediate result of part A and Theorem 4. Thus only part A 
needs to be proved. To proceed, note that |A+| < 1 and |A_| < 1 if and only if (i) 

(1-|A + | 2 ) (l — |A-| 2 ) > 0 

and (ii) 

(1-|A+| 2 |A_| 2 ) >0, 
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Thus, by using Eqs. (3.5), (2.4), (4.15), (4.20), and (4.21), it is easy to see that Eq. (4.2) 
is equivalent to Eq. (4.27) and 


r, s) — \ZeJp, t, s) >0, 0 < s < 1 


(4.28) 


As a result, to complete the proof, one needs only to show that Eqs. (4.25) and (4.26) is 
equivalent to Eq. (4.28). 

To proceed, for simplicity, in the following D(u,t,s ), E(u,t,s ), F(u, r, s), G(u,t,s), 
and s ) may be abbreviated as D, E, F, G, and H, respectively. By using the fact 

that E > 0 (see Eq. (4.19)), it is easy to show that the condition D — \f~E > 0 implies that 
(i) D > 0 and (ii) 

D 2 -E= (D + VE)(D-VE) >0 (4.29) 

Thus, with the aid of Eq. (4.22), one concludes that Eq. (4.28) implies both Eqs. (4.25) 
and (4.26). 

To show that Eqs. (4.25) and (4.26) imply Eq. (4.28), note that 


D - VE = D > 0 if T> > 0 and E = 0 (4.30) 

Moreover, because D + \[E > 0 if Zb > 0 and E > 0, one has 

jj2 jjj 

D — \J~E = — >0 if D > 0, D 2 — E > 0, and E > 0 (4.31) 

d + Ve ~ 


Thus, with the aid of Eqs. (4.19), (4.22), (4.30) and (4.31), one concludes that Eqs. (4.25) 
and (4.26) indeed imply Eq. (4.28). QED. 

At this juncture note that, given any (v. r), D(v, r, s), F(v, r, s) and G(v, r, s) are all 
quadratic polynomials in s and thus their minimum values in the interval 0 < s < 1 are 
easy to evaluate. As will be shown, this makes the analytical study of Eqs. (4.25)-(4.27) 
a relatively simple one. This is very fortunate because, according to Theorem 6, these 
equations play key roles in the current stability study. 

To proceed, note that an immediate result of Theorem 6 is Theorem 7. 

Theorem 7. (i) D(u,t, 0) > 0, (ii) D(v,t, 1) > 0, (iii) F(u,r, 0) > 0, (iv) F(v,t, 1) > 
0, (v) H(u,t, 0) > 0, and (vi) H(u,t, 1) > 0 are all necessary conditions for a given (u,t) 
to be c-t stable. 


To study conditions (i)-(vi) referred to above, Eqs. (4.16) (4.18), (4.23), and (4.24) 
are used to yield 

D(u,t, 0) = 4 (4.32) 

D(y, r, 1) = (2 - u 2 )r 2 + 2(2 - 3 u 2 )r + 2u 4 - bv 2 + 4 (4.33) 

F(u, r, 0) = 4 r (4.34) 
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and 


F{v, t, 1) = (2 + r + v 2 )(t - u 2 ) 
H(u, t, 0) = 0 


(4.35) 

(4.36) 


H(u, r, 1) = 4(1 - u 2 )(t - v 2 ) 2 [(2 + r) 2 - u 2 } (4.37) 

According to Eqs. (4.32) and (4.36), conditions (i) and (v) referred to in Theorem 7 are 
satisfied automatically. The significance of other conditions will be partially addressed in 
the following Theorems 8-11. 

Theorem 8. F(u, r, 0) > 0 and F(u. r, 1) > 0 if and only if r > v 2 . 

Proof. According to Eq. (4.34), F(u, r, 0) > 0 if and only if r > 0. With the aid of 
Eq. (4.35) and the fact that 2 + r + v 2 > 0 if r > 0, one concludes that F(u, r, 0) > 0 
and F(u, r, 1) > 0 imply r > v 2 . Conversely, it is easy to see that E(^, r, 0) > 0 and 
F( i/, r, 1) > 0 if r > v 2 . QED. 

Theorem 9. Let r > v 2 . Then H{v, r, 1) > 0 if and only if r > v 2 and v 2 < 1. 

Proof. With the aid of the assumption r > v 2 and Eq. (4.37), H{v , r, 1) > 0 implies 
(i) r > v 2 and (ii) 

( v 2 - 1) [u 2 - (2 + r) 2 ] > 0 (4.38) 

Because r > v 2 implies r > 0 and thus v 2 — 1 > v 2 — (2 + r) 2 , conditions (i) and (ii) 
imply either (a) v 2 < 1 or (b) v 2 > (2 + r) 2 . Case (b) can be ruled out because it along 
with condition (i) implies r > (2 + r) 2 , a result inconsistent with r > 0 which follows from 
condition (i). Thus H(u, r, 1) > 0 implies r > v 2 and v 2 < 1, if r > v 2 is assumed. 

Conversely, because (2 + r) 2 > r > v 2 if r > u 2 , Eq. (4.37) implies that H(u, r, 1) > 0 
if r > v 2 and v 2 < 1. Thus the proof is completed. QED. 

Theorem 10. Let r > v 2 . Then H{v , r, 1) = 0 if and only if at least one of the two 
cases: (i) r = v 2 and (ii) v 2 = 1, is true. 

Proof. Eq. (4.37) implies that H(u,r, 1) = 0 if and only if at least one of the three 
cases: (i) v 2 = 1, (ii) r = v 2 , and (iii) v 2 = (2 + r) 2 , is true. Case (iii) can be ruled out 
because it along with the assumption r > v 2 implies r > (2 + r) 2 , a result inconsistent 
with r > 0 (which follows from r >v 2 ). Thus the proof is completed. QED. 

Theorem 11. Let r = v 2 . Then D(u,r, 1) > 0 if and only if u 2 < 1. 

Proof. Let r = u 2 . Then Eq. (4.33) implies that 

D(u, t, 1) = (1 — t)(t 2 + 3r + 4) (r = u 2 ) (4.39) 

With the aid of Eq. (4.39) and the fact that 

r 2 + 3r + 4 = (r + 3/2) 2 + 7/4 > 7/4, — oo < r < +oo (4.40) 
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it is easy to see that, assuming r = u 2 , D(y , r, 1) > 0 if and only if v 2 < 1. QED. 

According to Theorems 8-10, the conditions (i) F{v,r, 0) > 0, (ii) F(u,t, 1) > 0, and 
(iii) H(u, r, 1) > 0 require that r = v 2 if the conditions r > v 2 and v 2 < 1 are not satished 
simultaneously. On the other hand, according to Theorem 11, the condition D(u, r, 1) > 0 
requires that v 2 < 1 for the case r = v 2 . Thus one has Theorem 12. 

Theorem 12. The conditions (i) D{ u, r, 1) > 0, (ii) F(u, r, 0) > 0, (iii) F(u , r, 1) > 0, 
and (iv) H(u,t , 1) > 0 require that t >v 2 and v 2 < 1. As such, Theorem 7 implies that 


r > v 2 and v 2 < 1 


are necessary conditions for a given {v. r) to be c-r stable. 


(4.41) 


In the following, it will be shown that only a subset of those r and v that satisfy 
the necessary conditions Eq. (4.41) will also satisfy the sufficient conditions for stability. 
As a prerequisite, we shall first study the conditions under which the matrix Q(is,t,6) is 
defective if r and v satisfy Eq. (4.41). We begin with Theorem 13. 


Theorem 13. Let r > v 2 and v 2 < 1. Then Q{v,t,6) is defective if and only if 


4(1 + r) = u 2 {t 2 + 2t + 5) 


and 


cos(0/2) = 0 


Proof. Assuming r > v 2 and v 2 < 1, first we will show that 

\ + (u,t, 9) = X-(u, t,6) 


(4.42) 

(4.43) 

(4.44) 


if and only if Eqs. (4.42) and (4.43) are satisfied. According to Eq. (4.10), Eq. (4.44) is 
equivalent to 

\[Y 2 + Y 2 + A = 0 and \[Y 2 + Y 2 - X = 0 (4.45) 


Thus Eq. (4.44) is true if and only if 


X = Y = 0 


(4.46) 


According to Eq. (4.6), Y = 0 if and only if at least one of the four cases: (a) v = 0, 
(b) r = 1, (c) sin(0/2) = 0, and (d) cos(#/2) = 0, is true. For case (a) v — 0, Eqs. (4.5) 
and the assumption r > v 2 imply that 


A = 4 [l + rsin 2 (#/2)] >4 
Thus case (a) is incompatible with Eq. (4.46). 


( u = o) 


(4.47) 
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For case (b) r = 1, Eq. (4.5) implies that 


X = 4cos 2 (#/2) + 8(1 — u 2 ) sin 2 (0/2) (r = 1) (4.48) 

Using the assumption i' 2 < 1, Eq. (4.48) implies that, for case (b), X = 0 if and only if 
v 2 = 1 and cos(#/2) = 0. 

Because cos 2 (#/2) = 1 if sin(0/2) = 0, Eq. (4.5) implies that X = 4 if sin(0/2) = 0. 
Thus case (c) is incompatible with Eq. (4.46). 

Because sin 2 (0/2) = 1 if cos(d/2) = 0, Eq. (4.5) implies that 

X = 4(1 + r) — v 2 {t 2 + 2t + 5) (cos(6>/2) = 0) (4.49) 

if cos(#/2) = 0. Thus, for case (d), X = 0 if and only if Eq. (4.42) is satisfied. 

Assuming r > v 2 and v 2 < 1, it has been shown that X = Y = 0 if and only if at 
least one of the following two conditions: (i) 

r = 1, u 2 = 1, and cos(d/2) = 0 (i.e., case (b)) 


and (ii) 


cos(#/2) = 0 and 4(1 + r) = u 2 (r 2 + 2r + 5) (i.e., case (d)) 

is met. Because r = 1 and v 2 = 1 form a special solution of Eq. (4.42), condition (i) is 
only a special case of condition (ii). Thus, assuming r >i/ 2 and J' 2 < 1, Eq. (4.44) (which 
is equivalent to X = Y = 0) is true if and only if Eqs. (4.42) and (4.43) are satisfied. 
Moreover, with the aid of Eq. (2.8) and the fact that sin(0/2) = ±1 if cos(#/2) = 0, 
Eq. (4.43) also implies that one of the off-diagonal elements of Q(u, r, 9) does not vanish 
and thus Q(y , r, 9) is not a multiple of /. According to Theorem 1, Q(y , r, 9) is defective if 
and only if (i) Eq. (4.44) is true and (ii) Q(v, r, 9) is not a multiple of /. Thus the current 
theorem is proved. QED. 

An immediate result of Theorem 13 is Theorem 14. 

Theorem 14. The matrix Q(u, r, 9) is defective if r = v 2 = 1 and cos(#/2) = 0. 

To proceed, we will establish Theorem 15. 

Theorem 15. Let Q{y , r, 9) be defective with r > v 2 and v 2 < 1. Then the special 

case 

p(Q(v,r,9)) = 1 (4.50) 

occurs if and only if 

r = v 2 = 1, and cos(#/2) = 0 (4-51) 

Proof. As a preliminary, first we will deduce several results from the current basic 
assumption, i.e., Q{y,r,9) is defective with r > and < 1. According to Theorem 13 
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and its proof, Eqs. (4.42), (4.43), and (4.46) follow immediately from the basic assumption. 
Also, by using Eq. (4.42) and the fact that 

t 2 + 2t + 5 = (1 + t ) 2 + 4 > 4, — oo < r < +oo (4-52) 


one concludes that 

2 = 4 (1 + r ) 

r 2 + 2r + 5 


(4.53) 


Moreover, because sin(0/2) 
Eq. (4.10) implies that 


±1 if cos(d/2) = 0, with the aid of Eqs. (4.43) and (4.46), 


p(Q(u, t, 9)) 


u(S + r) 

2(1+7) 


(4.54) 


Next assume Eq. (4.50). Because 3 + r > 0 (which follows from the assumption 
r > i/ 2 ), Eqs. (4.50) and (4.54) imply that 


.2 4(1 +T) 2 

(3 + t) 2 


(4.55) 


Eliminating ;/ 2 from Eqs. (4.53) and (4.55) and using the basic assumption Eq. (2.4) (which 
is consistent with the current assumption r > u 2 ), one has 


r 3 + 2 r 2 + r — 4 = (r — l)(r 2 + 3r + 4) = 0 


(4.56) 


Eq. (4.56) coupled with Eq. (4.40) implies that r = 1. In turn, by using either Eq. (4.53) 
or Eq. (4.55), one has v 2 = 1 as a result of r = 1. Because Eq. (4.43) (i.e., cos(#/2) = 0) is 
a result of the basic assumption, it has been shown that Eq. (4.51) follows from the basic 
assumption and Eq. (4.50). 

Conversely, with the aid of (i) Theorem 1, and (ii) Eqs. (2.8) and (3.5), it can be 
shown by direct substitution that both the basic assumption and Eq. (4.50) are valid for 
the special case Eq. (4.51). Thus the proof is completed. QED. 

Next we have Theorem 16. 

Theorem 16 . A given ( v , r) satisfies Eq. (4.2) and yet is c-r unstable if and only if 
r = v 2 = 1. 

Proof. Theorems 6 and 12 imply that Eq. (4.41) is a result of Eq. (4.2). Thus, 
according to Theorems 5 and 15, r = v 2 = 1 if {v. r) satisfies Eq. (4.2) and is also c-r 
unstable. 

Conversely, Theorem 6 coupled with Eqs. (4.16), (4.18), and (4.23) implies that any 
(u, t) with r = v 2 = 1 satisfies Eq. (4.2). Moreover, according to Theorems 3, 14 and 15, 
such a ( v , r) is also c-r unstable. Thus the proof is completed. QED. 

At this juncture, note that Theorems 14 and 15 state that, for the special case 
Eq. (4.51), Q{y,r,6) is defective with p(Q(u,r,9)) = 1. Thus, according to a com- 
ment made following Eq. (3.16), for this special case, the magnitude of any element in 
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[Q{y, r, 6)] m will grow not faster than linearly with to. Because round-off errors associated 
with a modern computer are in the order of 10 -10 or less, the instability associated with 
this special case generally is very mild and may not be detected even after billions of time 
steps have elapsed. 

Next, by combining Theorems 6, 12 and 16, one arrives at Theorem 17. 

Theorem 17. A given (v, r) which does not satisfy Eq. (4.41) is c-r unstable. On 
the other hand, a given (u, r) which satisfies Eq. (4.41) is c-r stable if and only if (i) it 
satisfies Eqs. (4.25)-(4.27); and (ii) it does not belong to the special case r = u 2 = 1. 

Compared to those given in Theorem 3, the necessary and sufficient stability conditions 
given in Theorem 17 are much more explicit and easier to handle. As such, this theorem 
will be used repeatedly in the rest of the development. In particular, it will be used to 
establish Theorem 18. 

Theorem 18. The c-r scheme is stable for any one of the following special cases: (a) 

v = 0 and r > 0; (b) v 2 = 1 and r > 1; and (c) 0 < v 2 < 1 and r = \v\. 

Proof. Let 0 < s < 1 throughout this proof. Then, with the aid of Eqs. (4.16), (4.18), 
(4.23), and (4.24), for case (a) u = 0 and r > 0, one has 

D(v , r, s) = D( 0, r, s) = 2 [(1 + rs) 2 + l] > 4 (4-57) 

F( v, r, s) = F{ 0, r, s) = r(2 — s)(2 + rs) > 0 (4.58) 

and 

H{y , r, s ) = H( 0, r, s) = 4s 2 r 2 (2 + rs) 2 > 0 (4.59) 

Because v = ±1 if v 2 — 1, for case (b) v 2 = 1 and r > 1, one has 

D(u, r, s) = D(± 1, r, s) = (1 — r) 2 s + 4(1 — s) > 0 (4.60) 

F(y, r, s) = F(±l, r, s) = (1 — r) 2 s + 4(r — s) > 0 (4.61) 

and 

H(u,t,s) = H(±1,t,s) = 0 (4.62) 

Because 0 < v 2 < 1 and r = \u\ if and only if v = ±r and 0 < r < 1, for case (c) 
0 < v 2 < 1 and r = \v\, one has 

D(v, r, s) = D(±t, r, s) = r(l — r)(8 + 5r — r 2 )s + 4(1 — rs) > 0 (4.63) 

F(u, r, s) = E(±r, r, s) = r(l — r)(r 2 + r + 2)s + 4r(l — s) > 0 (4.64) 

and 

H{v, r, s) = H(±t, r, s) = 16r 2 (l — r 2 ) 2 (l — r)s 2 > 0 (4.65) 
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Obviously cases (a) and (b) are special cases of the more general case defined by 
Eq. (4.41). Moreover, because v 2 < \v\ if 0 < v 2 < 1, case (c) is also a special case of the 
more general case. In addition, none of cases (a)-(c) contains the special case r = v 2 = 1. 
With the aid of these observations and Eqs. (4.57)-(4.65), Theorem 18 follows directly 
from Theorem 17. QED. 


Next let 


fit = {(zq t)|0 < v 2 < 1,t > v 2 and r 2 ^ z, 2 } (4.66) 

4'_ = f {(zq r)|0 < v 2 < 1,t > v 2 and r 2 < z^ 2 } (4.67) 

and 

T + = f {(u, r)|0 < v 2 < l,r > v 2 and r 2 > u 2 } (4.68) 

Then and are disjoint, and 

$ = E + UE_ (4.69) 


Moreover, we have Theorems 19 and 20. 

Theorem 19. Excluding the four special cases addressed in Theorems 16 and 18, 
T is the set of all other (v. r) that satisfy the necessary stability conditions r > v 2 and 
v 2 < 1 given in Theorem 12. 

Proof. Note that (i) r = \v\ > v 2 if 0 < z' 2 < 1 and r = \v\\ (ii) r 2 = z' 2 if r = |z/|, 
(iii) t = \p\ ii t > p 2 and r 2 = u 2 , and (iv) r = r 2 = u 2 implies either r 2 = u 2 = 0 
or t 2 = v 2 = 1. Items (i)-(iii) imply that 0 < z^ 2 < 1 and r = \u\ (which is case (c) in 
Theorem 18) if and only if 0 < z^ 2 < 1, r > z^ 2 , and r 2 = v 2 . On the other hand, item (iv) 
implies that the case with both 0 < v 2 < 1 and r = r 2 = v 2 does not exist. The proof 
follows from the above two observations and the facts that (i) r > v 2 = 0 if and only if 
v = 0 and r > 0, and (ii) t > u 2 = 1 ii and only if either (a) r = v 2 = 1 or (b) v 2 = 1 and 
r > 1. QED. 

Theorem 20. Eq. (4.68) is equivalent to 

\H_i_ = {(z^, r)|0 < v 2 < 1, r > v 2 and r 2 > z' 2 } (4.70) 

Proof. Note that (i) u 4 > v 2 if r = v 2 and r 2 > u 2 , and (ii) the relations u 4 > v 2 and 
0 < v 2 < 1 are contradictory. Thus the case with 0 < v 2 < 1, r = u 2 , and r 2 > v 2 does 
not exist, i.e., Eq. (4.68) is equivalent to Eq. (4.70). QED. 

To proceed, we will establish Theorems 21 and 22. 

Theorem 21. Let (zq t) G fib Then 

D(p, r,s) > 0, 0 < s < 1 (4-71) 
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Proof. As a preliminary, note that Eq. (4.33) implies that 


D(u,t, 1) = (2- v 2 ) 



2-3 z/ 2 \ 

2-i/ 2 ) 


2(1 — z/ 2 )(z/ 4 + v 2 + 2) 
(2 — 1 / 2)2 


z/ 2 ^ 2 (4.72) 


Thus 

,D(z/, r, 1) > 0 if i/ 2 < 1 (4.73) 

Let (z/, r) G T_. Then Eqs. (4.16) and (4.67) imply that 

a 2 g( y ,r.5) l =4(1 _^ 2 )(t 2 _ i/ 2 )<0 ((i 2 ,r) € *&-) (4.74) 

J I/,T 

i.e., for any given (z/, r) € T_, the relation between the function D(u, r, s) and s is repre- 
sented by a curve which is concave downward on the s-D plane. Thus 

min D(v, r, s) = mm{D(v, r, 0), D(u, r, 1)} ((z/, r) 6 $_) (4-75) 

0<s<l 

By using Eqs. (4.32) and (4.73), Eq. (4.75) implies that 

D(is,t,s)> 0, 0<s<l ((z/, r) G T_) (4.76) 

Next let (u, r) G T + . Then, by using Eq. (4.68) (in particular the facts that z/ 2 < 1 
and (1 — z/ 2 )(r 2 — u 2 ) > 0), Eq. (4.16) implies that 


D(v, r, s) > [4t + (r 2 — 6 r — 3) z/ 2 ] s + 4 > [4 rz/ 2 + (r 2 — 6r — 3)z/ 2 ] s + 4 
= (1 — r) 2 z/ 2 s + 4(1 — z/ 2 s) >0, 0 < s < 1 ((z/, r) G T + ) 


(4.77) 


It has been shown that D(u,r,s ) > 0, 0 < s < 1, for both case (a) (z/, r) G and 
case (b) (u, r) G T + . Because T = U T + , the proof is completed. QED. 

Theorem 22. Let (z/, r) G T. Then 


F(u, t,s) > 0, 0 < s < 1 


(4.78) 


Proof. Let (u,t) G T + . Then Eqs. (4.18) and (4.68) imply that 

a 2 F(^T,s) l =2(1 _^2 )(;/ 2_ r 2 )<0 ((i/.rJe'M (4.79) 

J V,T 

i.e., for any given (z/, r) G the relation between the function F(u,t,s ) and s is repre- 
sented by a curve which is concave downward on the s-F plane. Thus 

min F(u, r, s) = mm{F(u, r, 0), F(u, r, 1)} ((z/, r) G f+) (4.80) 


NASA/TM— 2005-2 1 3627 


20 



By using Eqs. (4.34), (4.35) and (4.70), Eq. (4.80) implies that 

F(is,t,s)> 0, 0<s<l ((is, r) G \l/_)_) (4-81) 

Next let ( v,t ) G T_. Then, by using Eq. (4.67) (in particular the facts that (1 — 
v 2 )(is 2 — t 2 ) > 0 and 0 < r < \v\ < 1), Eq. (4.18) implies that 


= 2(1 - is 2 )(is 2 - t 2 )s - [2r(l - r) + (3 + t 2 )v 2 ] 

V,T 

< 2(1 - is 2 )(v 2 - T 2 ) - [2r(l - r) + (3 + t 2 )p 2 ] (4.82) 

= -2(1 - is 2 )t 2 - 2is 4 - 2r(l - r) - (1 + t 2 )is 2 < 0, 

0 < s < 1 ((is, t) G T_) 

Thus, for any given (u,t) G T_, the relation between F and s is represented by a curve 
on the s-F plane which has a negative slope in the interval 0 < s < 1. In turn, this fact 
coupled with Eqs. (4.35) and (4.67) implies that 

F(is, t, s ) > F(is, r, 1) > 0, 0 < s < 1 ((is, r) G T_) (4.83) 

It has been shown that F(v,r,s) > 0, 0 < s < 1, for both case (a) (v,t) G and 
case (b) (is,r) G T_. Because T = U T + , the proof is completed. QED. 

According to Theorems 21 and 22, Eqs. (4.25) and (4.27) are satisfied by all (is, r) G T. 
Thus, Theorem 17 implies that a given (is, r) G T is c-r stable if and only if it satisfies 
Eq. (4.26). Thus, with the aid of Eqs. (4.23) and (4.66), one arrives at Theorem 23. 

Theorem 23. For any given (is,t) G T, Eq. (4.26) is equivalent to 

inf G(is, t,s) > 0 (4.84) 

0<s<l 


dF(v, r, s) 

ds 


where the expression on the left side of the sign “>” denotes the inhmum (i.e., the greatest 
lower bound) of G(is,r,s ) in the interval 0 < s < 1. As such, a given (is,t) G T is c-r 
stable if and only if it satisfies Eq. (4.84). 

Because of Theorem 23, in the following we shall focus on finding those (is, r) G T that 
satisfy Eq. (4.84). 

To proceed, first we will establish Theorem 24. 

Theorem 24. For any given (is,r) G T, let 


s Q (is, t) 


def u 2 t 2 + (4 - 6is 2 )t - 3 v 2 

~ 2(1 — z/ 2 )(iA — t 2 ) 


(4.85) 
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Let s 0 (u,t) be abbreviated as s Q . Then 


inf G(u,t,s) 

0<s<1 


G(u, r, s a ) if 0 < s Q < 1 
G(u, r, 1) if s 0 > 1 
G(u, r, 0) if s Q < 0 


(4.86) 


Proof. To facilitate the proof, the domain of the function G defined in Eq. (4.24) 
will be extended to — oo < s < +oo. As such, for any given (is, r) G T and any s with 
— oo < s < Too, one has 


dG(u, t, s ) 

ds 


V,T 


2(1 - v 2 ){r 2 - v 2 ) 2 [s - s Q (u, r)] 


(4.87) 


and 

= 2(1 - v 2 ){t 2 - v 2 ) 2 > 0 (4.88) 

Thus, for any given (u,t) G T. (i) the relation between the function G(u,r,s ) and s is 
represented by a curve which is concave upward on the s-G plane, and thus the absolute 
minimum of G in the interval — oo < s < Too occurs at where dG/ds = 0, i.e., 


d z G(u 1 r, s) 
d? 


s = s 0 (u, t) 


(4.89) 


(ii) G is strictly monotonically decreasing in the interval s < 1 if s G > 1: and (iii) G is 
strictly monotonically increasing in the interval s > 0 if s Q < 0. In addition, for any given 
(zy t), because G is a continuous function of s in the interval — oo < s < Too, one also has 

( iv ) 

lim G(u,t, s) = G(u,t, 0) (4.90) 

S^0 + 

Eq. (4.86) is a direct result of (i)-(vi). QED. 

With the aid of Theorem 24, the bulk of the remaider of the paper will be devoted to 
answer a key question, i.e., given any v with 0 < v 2 < 1 (which is required by the condition 
(u, t) G dl), what is the range of r that will satisfy Eq. (4.84) and the rest of the condition 
(u, r) G \l/ (i.e., t>v 2 and r 2 ^ v 2 )l 
To proceed, let 


I±(x) 


def 3a: — 2 T 2y/3x z -3x + l 

x 


and (iii) 


T / \ def 3a: 

J ± [x) = — 


2 T ^2(a: 3 
2 — x 


x + 2) 


0 < x < 1 


0 < x < 1 


(4.91) 


(4.92) 
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Hereafter, for any function f(x), as usual \/f(x) denotes the principal square root of f(x). 
As such f{x) > 0 if f(x) > 0. Given Eqs. (4.91) and (4.92), one can establish Theorem 
25. 

Theorem 25. In the domain 0 < x < 1, we have 



I+(x ) > 0 

(0 < x < 1) 

(4.93) 


I-(x ) < 0 

o 

A 

A 

(4.94) 


J+(x ) > 0 

(0 < X < 1) 

(4.95) 

and 

J- ( x ) < 0 

(0 < x < 1) 

(4.96) 

Proof. Because 

A(3x 2 - 3x + 1) = 

(3x - 2) 2 + 3^ 2 

(4.97) 

one has 

2\j3x 2 — 3x + 1 > \3x — 2 1 , x ^ 0 

(4.98) 

Eqs. (4.93) and (4.94) follow directly from Eqs. (4.91) and (4.98). 
Next because 



to 

CO 

1 

+ 

II 

1 

! + 2x(x -2) ^x - ^ 

(4.99) 

one has 

\/2(x 3 — x + 2) > |3x 

— 2|, 0 < x < 2 

(4.100) 


Eqs. (4.95) and (4.96) follow directly from Eqs. (4.92) and (4.100). QED. 

With the above preparations and the understanding that hereafter the symbol “<G>” 
may be used to take the place of the statement “if and only if” , Theorem 26 can now be 
presented. 

r) G T_, we have 
5=> r > I + {v 2 ) 

t = I + (u 2 ) ((^,r)GT_) (4.101) 

5=> r < I + {v 2 ) 

^ r > J+(v 2 ) 

^ r = J + (u 2 ) (4.102) 

r < J + {y 2 ) 


Theorem 26. (A) For any (u, 

> 0 



and 
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On the other hand, (B) for any (is, r) G \b_|_, we have 


and 


s 0 (z/, r) 


> 0 

<^> 

r < I+(is 2 ) 

= 0 


r = I + (u 2 ) 

< 0 


T > I + (lS 2 ) 

> 1 


T < J+(lS 2 ) 

= 1 

<^> 

II 

£ 

to 

< 1 


r > J + (is 2 ) 


((z/,r) e tf+) 


((j/,t)g*+) 


(4.103) 


(4.104) 


Proof. As a preliminary, note that 

u 2 r 2 + (4 — 6p 2 )t — 3v 2 = v 2 [t — I + (u 2 )] [r — I-(is 2 )] (0 < v 2 < 1) 

In addition, because r > v 2 and 0 < v 2 < 1 if (is, r) G \b, Eq. (4.94) implies that 

T — I-{y 2 )> 0, (zq r) G 4/ 


(4.105) 


(4.106) 


Because the expression on the left side of Eq. (4.105) is the numerator of the fraction 
on the right side of Eq. (4.85), Eq. (4.101) now follows from Eqs. (4.85), (4.105) and 
(4.106), and the fact that 0 < v 2 < 1, and v 2 — r 2 > 0 if (is, t) G 

To prove Eq. (4.102), note that Eq. (4.85) implies that, for any (is, r) G 4/, 


s 0 (zqr) - 1 


(2 — z/ 2 )r 2 + (4 — 6z A )t — zA(5 — 2zA) 
2(1 — z/ 2 )(z/ 2 — r 2 ) 


(4.107) 


Also one has 

(2 - z/ 2 )r 2 + (4 - 6is 2 )t - z/ 2 (5 - 2z/ 2 ) = (2 - z/ 2 ) [r - J + (z/ 2 )] [r - J-(is 2 )] (0 < z/ 2 < 1) 

(4.108) 

In addition, because r > z/ 2 and 0 < z/ 2 < 1 if (is, r) G 41, Eq. (4.96) implies that 

r — J-(is 2 ) > 0, (is, r) G lb (4.109) 

Because the expression on the left side of Eq. (4.108) is the numerator of the fraction 
on the right side of Eq. (4.107), Eq. (4.102) now follows from Eqs. (4.107)-(4.109), and 
the fact that 0 < v 2 < 1 and v 2 — r 2 > 0 if (is, r) G lb-. 

This finishes the proof of part A. Part B can be proved using a fine of logic identical 
to that used to prove part A. The only difference that sets part B apart from part A is 
that is 2 — t 2 < 0 for the case (is, r) G 4/ + while is 2 — r 2 >0 for the case (is, r) G 4/_ . QED. 
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Next, note that Eq. (4.24) yields 


G(u, r, 1) = (r — u 2 ) 2 [(2 + r) 2 — v 2 ] 


and 

G(u, r, 0) = 4 r [u 2 t 2 + (1 - v 2 )t - v 2 ] 
In addition, for any (u,t) £ 'll, Eqs. (4.24) and (4.85) also yield 

,2 


G(v, t, s 0 ) = 


u“(l + t) 2 [u 2 t 2 + 2 (y 2 — 4 )t + 9u 2 ] 


4(1 - v 2 ) 

An immediate result of Eqs. (4.66) and (4.110) is Theorem 27. 
Theorem 27. For any (u, r) £ T, we have 

G(u,t, 1) > 0 ((i/,r)eT) 


Next let 


K ± (x) 


def x — 1 4 Vl — 2x + 5x 2 
2x 


0 < x < 1 


Then one has Theorems 28 and 29. 

Theorem 28. In the domain 0 < x < 1, we have 


and 

Proof. Because 
one has 


K + (x) >0 (0 < x < 1) 

K-(x) <0 (0 < x < 1) 

1 — 2x + hx 2 = (x — l) 2 + 4a; 2 


yjl — 2a; + 5a; 2 > |ar — 1 1 , x^0 

Eqs. (4.115) and (4.116) follow directly from Eqs. (4.114) and (4.118). QED. 
Theorem 29. For any (u,r) £ T, we have 

G(u,t, 0) > 0 r > K + (u 2 ) ((u,r)£'ff) 

Proof. Note that 

4 r [u 2 t 2 + (1 - p 2 )t - v 2 ] = 4 tv 2 [t - K + (u 2 )] [r - K_(v 2 )\ , 0 < v 2 < 1 


(4.110) 

(4.111) 

(4.112) 

(4.113) 

(4.114) 

(4.115) 

(4.116) 

(4.117) 

(4.118) 

(4.119) 

(4.120) 
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In addition, because r > v 2 and 0 < v 2 < 1 if (z/, r) G \H, Eq. (4.116) implies that 

r — K-(v 2 ) > 0, (4.121) 

Eq. (4.119) now follows from Eqs. (4.111), (4.120) and (4.121), and the fact that r > v 2 
and 0 < v 2 < 1 if (v. r) G 4/. QED. 

Next let 


r , , def 4 — X ± 2y/2(2 — X — X 2 ) 
T±(x) = 


0 < x < 1 


X 


(4.122) 


Then one has Theorems 30 and 31. 

Theorem 30. In the domain 0 < x < 1, we have 


L + (x) > L-(x) > 0 (0 < x < 1) (4.123) 


Proof. Note that (i) 

2 — x — x 2 = — (x + 2)(x — 1) > 0, —2 < x < 1 


and (ii) 
Thus 


(4 — x) 2 — 2y/2(2 — x — x 2 ) = 9x 2 > 0, x^0 


4 — x = |4 — x| > 2y/2(2 — x — x 2 ) >0, 0<x<lor— 2<x<0 

Eq. (4.123) is a result of Eqs. (4.122) and (4.126). QED. 

Theorem 31. For any (u,t) G dl, we have 

G(v,t,s 0 )> 0 L_{v 2 ) < r < L + (is 2 ) ((i/, r) G W) 


(4.124) 


(4.125) 


(4.126) 


(4.127) 


Proof. Note that 

z/ 2 r 2 + 2(z/ 2 - 4)r + 9i/ 2 =z/ 2 [r-L + (i/ 2 )] [r-L_(z/ 2 )] (0 < zv 2 < 1) (4.128) 

Because 1 + r > 0, r/ 2 > 0, and 1 — z^ 2 > 0 if (v,t) G dl, Eqs. (4.112) and (4.128) imply 
that 

G{v, r, s 0 ) > 0 [r - L + (z/ 2 )] [r - L_(z/ 2 )] < 0 ((zqr)GT) (4.129) 


NASA/TM— 2005-2 1 3627 


26 



if (v,t) G T. Because 0 < v 2 < 1 if (u,r) G 4/, Eq. (4.127) now follows from Eq. (4.129) 
and a result of Eq. (4.123), i.e., 

[t — L_(i/ 2 )] > [t — L + (i/ 2 )] , 0 < z/ 2 < 1 (4.130) 

QED. 

With the above preliminaries, one can establish Theorem 32. 

Theorem 32. (A) Let (u, r) G Then (u, r) is c-r stable if and only if it satisfies 
one of the three mutually exclusive sets of conditions specified, respectively, in Eqs. (4.131)- 
(4.133): 

r > J+(u 2 ) (4.131) 

K+(u 2 ) < r < I + (v 2 ) (4.132) 

and 

I+(u 2 ) <r < J+(p 2 ) and L-(v 2 ) < t < L + (v 2 ) (4.133) 

(B) Let (v,t) G T_|_. Then [u. r) is c-r stable if and only if it satisfies one of the three 
mutually exclusive sets of conditions specified, respectively, in Eqs. (4.134)-(4.136): 

r < J+(u 2 ) (4.134) 

r > I+{v 2 ) and r > K + (is 2 ) (4.135) 

and 

J_i_(z/ 2 ) < r < /+(^ 2 ) and L_(i/ 2 ) < r < I/ + (i/ 2 ) (4.136) 

Proof. Let 

^(a) def ,r^ T )|^ T ) g v[f_ an( J So (v,t) > 1} (4.137) 

^(d) d| f G and s 0 (u,t) < 0} (4.138) 

^( 7 ) d£ f |(^ )T )|(^ )T ) ^ 4/_ and 0 < s 0 (j/, r) < 1} (4.139) 

'L+) = f {(u,t)\(u,t) G and s 0 (v,t) > 1} (4.140) 

^( d ) d£ f {( j , } r )|( j , }r ) ^ \ j / + and s 0 (u,t) < 0} (4.141) 

and 

4/(7) = f {(u,t)\(u,t) G \b_i- and 0 < s 0 (u,r) < 1} (4.142) 

Because \P_ and T_|_ are mutually exclusive, the above definitions imply that (i) 4/^, 
^,(d), xj/L) , , and are mutually exclusive; (ii) 

T_ = ^ U U (4.143) 
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and (iii) 



\& + = yW u U ^ 

(4.144) 

Moreover, by using Theorem 26, Eqs. (4.137)-(4.142) imply 



_ {(v,t) \(is,t) E 4/- and r > J+(is 2 )} 

(4.145) 


^(d) _ t)\(is, t) E 41 _ and r < J_|_(^ 2 )} 

(4.146) 


4/ ( _ 7) = {(is,r)\(is,r) E and I+(v 2 ) < r < J+(is 2 )} 

(4.147) 


= {(u,t)\(u,t) E 4'+ and r < J+(u 2 )} 

(4.148) 


= {(v,t)\(v,t) E W + and r > I+(v 2 )} 

(4.149) 

and 

= {(is,t)\(is,t) E 4^ + and J+(is 2 ) < r < I+(u 2 )} 

(4.150) 


respectively. 

To proceed, note that: 

(a) With the aid of (i) Eqs. (4.137) and (4.140), and (ii) Theorems 24 and 27, Theorem 
23 implies that a given (is, r) E 4/^ U 4/^ is always c-r stable. 

(b) With the aid of (i) Eqs. (4.138) and (4.141), and (ii) Theorems 24 and 29, Theorem 
23 implies that a given (is, r) E 4'^^ U 4 /l ^ ) is c-r stable if and only if 

t>K + (u 2 ) (4.151) 

(c) With the aid of (i) Eqs. (4.139) and (4.142), and (ii) Theorems 24 and 31, Theorem 
23 implies that a given (is, r) E 41^ U is c-r stable if and only if 

L-(u 2 ) <r <L+(v 2 ) (4.152) 

Theorem 32 now follows from Eqs. (4.143)-(4.150) and the facts presented in the above 
items (a)-(c). QED. 

In principle, the question of whether a given (is, r) is c-r stable can now be answered 
by using Theorems 12, 16, 18, 19, and 32. However, in its current complicated form, 
Theorem 32 is difficult to use. Fortunately, Theorem 32 can be simplified greatly and, in 
fact, the stability condition for the c-r scheme can be cast into a rather simple explicit 
form. To obtain this simple form, we begin with Theorem 33. 

Theorem 33. We have: (A) 

(is, r) € 4/_ 0 < is 2 < 1 and is 2 < r < (4.153) 
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(B) 'I'- is not empty; and (C) 


(u, r) G T + 0 < v 2 < 1 and r > Vv 2 


(4.154) 


Proof. Because (i) —Vu 2 < r < \/~J 2 if r 2 < i/ 2 , and (ii) r 2 < v 2 if 0 < r < Vv 2 , part 
A is an immediate result of Eq. (4.67). Part B follows from the trivial fact that v 2 < Vu 2 
if 0 < v 2 < 1. To prove part C, note that (i) r > 0 if v 2 > 0 and r > u 2 , and (ii) r > Vv 2 
if r > 0 and r 2 > v 2 . Thus Eq. (4.70) implies that 0 < z^ 2 < 1 and r > Vu 2 if (u, r) G T+. 
Conversely, because (i) v^ 2 > v 2 if 0 < v 2 < 1; (ii) r > v 2 if r > 'Jv 2 and Vv 2 > v 2 \ 
and (iii) r 2 > v 2 if r > V^ 2 , one concludes that (u, r) G if 0 < v 2 < 1 and r > Vv 2 . 
QED. 

Next let 


def „ 
Ci = 3 


2\/2 


c 2 d = 3/11 


d(3f / . 1 

c 3 = (41 


7y / 33)/2 


and 


def 
C 4 = 


1664 181 

+ 


1664 181 


27 


27 


27 


27 


We have (i) c 4 ~ 0.172, c 2 ~ 0.273, C3 ~ 0.394 and c 4 ~ 0.530, and (ii) 

0 < ci < c 2 < c 3 < c 4 < 1 


(4.155) 

(4.156) 

(4.157) 

(4.158) 


(4.159) 


With the above preparations, we have Theorem 34. 

Theorem 34. (A) In the domain 0 < x < 1, / + (x) . J+ (x), K + (x), and L-(x) are 
strictly monotonically increasing while L + (x) is strictly monotonically decreasing; (B) we 
have 


I+(x) < x < K + (x) < L-{x) < J+(x) < < L + (x), 

I+(x) = x < K + (x ) < L_(x) < J+(x) < \fx < L + (x), 
x < I+(x) < K+(x) < L_(x) < J+(x) < \fx < L + (x), 
x < I+(x) = K + (x ) = L-(x) < J+(x) < \fx < L+(x), 
x < K + (x) < L_(x) < I+(x) < J+(x) < \fx < L + (x), 
x < K + (x) < L-(x) < I+(x) = J+(x) = yfx < L+(x), 
x < K+(x) < L_(x) < \fx < J+(x) < I+(x) < L + (x), 


0 < x < ci 

(4.160) 

x = c 1 

(4.161) 

Cl < X < c 2 

(4.162) 

X = c 2 

(4.163) 

C 2 < X < c 3 

(4.164) 

X = C2, 

(4.165) 

C 3 < X < c 4 

(4.166) 
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x < K + (x) < L-(x) = \fx < J+{x) < I+(x) < L+(x), x = C 4 (4.167) 

and 

x < K + (x) < \[x < L-(x) < J+(x ) < I+(x ) < L + (x), C4 < x < 1 (4.168) 

(C) 

K' + (c 2 ) = L'_(c 2 ) = 121/90 (4.169) 

where K' + (x) = f dK + (x) /dx and L'_(x) dL-(x)/dx] and (D) 

lim L_(x) = 0 and lim K + (x) = 1 (4.170) 

#^ 0 + x^l~ 

In order not to interrupt the current stream of development, the lengthy proof for 
Theorem 34 will be provided later in the paper. Here, with the aid of this theorem, we 
shall establish a simplified form of the stability condition for the c-r scheme as given in 
Theorem 35. 

Theorem 35. Let 


7 0 if x = 0 


T 0 (x) 


def 


L-{x) if 0 < £ < 3/11 
K + (x) if 3/11 < x < 1 


l 1 if x = 1 


(4.171) 


To = {(u,t)\u 2 < 1 , T > T 0 {y 2 ) and ( u 2 ,r ) ^ (1, 1)} (4.172) 


T = {(u,r)\u 2 < 1 and r > r 0 (u 2 )} (4.173) 

Then: (A) t 0 (x) is continuous at x = 0 and x = 1; (B) r 0 (x) is consistently defined at 
x = 3/11; (C) 

lim r'(x) = lim t'Jx) = 121/90 (4.174) 

™ . _ 3 _ — _ , 3 + 


where r'(x) f = dr 0 (x) / dx: (D) r 0 (x) is strictly monotonically increasing in the interval 
0 < x < 1; (E) 

x < t 0 (x) < \fx. 0 < x < 1 (4.175) 

(F) a given (v, r) is c-r stable if and only if (u, r) € T 0 ; and (G) a given (v, r) satisfies 
Eq. (4.2) if and only if (v. r) e T. 

Proof. Part A is a result of Eqs. (4.170) and (4.171). Part B follows from the fact 
that L_(3/ll) = K + (S/11) = 1/3. Part C follows from Eqs. (4.156) and (4.169). Part D 
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is a result of part A of Theorem 34, and parts B and C of the current theorem. Part E is 
a result of Eqs. (4.160)-(4.168) and (4.171). 

To prove part F, one needs to show that: (i) (u, r) G T 0 for any (u, r) that is c-r 
stable; and (ii) (u, r) ^ T 0 for any (v. r) that is c-r unstable. Here whether any particular 
(zy r) is c-r stable is determined using Theorems 12, 16, 18, 19, and 35. 


To proceed, let 



‘= {(z/, T )\v 2 > 1 or r < v 2 < 1 } 

(4.176) 


$2 = f {{v,t)\t = v 2 = 1 } 

(4.177) 


$3 *== {(v,t)\v 2 = 0 and r > 0 } 

(4.178) 


$4 *== {( v,t)\v 2 = 1 and r > 1 } 

(4.179) 


$5 d = {(z', r) 0 < v 2 < 1 and r = \u\} 

(4.180) 

With the aid Theorem 19, it is seen that T_, and the five sets defined above are 

inclusive and yet mutually exclusive, i.e., any (zy r) belongs to one and only one of these 
sets. To facilitate the proof, and T+ , respectively, will be further divided into several 

disjoint subsets to be defined immediately. 

Let 

^,( 1 ) A) t)|0 < v 2 < C 2 and v 2 < r < \/zA} (4.181) 


^(2) def ^2 _ C2 an( j j/2 < t < \/zA} 

(4.182) 


^>(3) <l| f C2 < z / 2 < C 3 and v 2 <r < 

(4.183) 

and 

^,(4) r)|c 3 < v 2 < 1 and u 2 <r < VzA} 

(4.184) 

Because (z/, r) G 4/ 
that (i) 4'^, £ = 1 

_ 0 < z / 2 < 1 and v 2 < r < \/zA (see Theorem 33), one 

, 2, 3, 4, are nonempty disjoint subsets of T_, and (ii) 

concludes 


= u} =1 *<? 

(4.185) 

Next let 

^(i) def 'y— ) 1 0 < v 2 < C 3 and r > •s/z' 2 } 

(4.186) 

and 

^( 2 ) def r)|c 3 < is 2 < 1 and r > \/zA} 

(4.187) 

Because (u, r) G \E. 
and are 

f 4^* 0 < v 2 < 1 and r > \/zA (see Theorem 33), one concludes that (i) 
nonempty disjoint subsets of T_)_, and (ii) 


W + = 4^ U T^ 2) 

(4.188) 
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From the above discussion, the sets (i) 4q, £ = 1, 2, 3, 4, 5; (ii) £ = 1,2, 3, 4; and 

(iii) 4^ and 4^ , are inclusive and yet mutually exclusive, i.e., any (za r) must belong to 
one and only one of these sets. Part F will be proved by showing that it is valid over each 
of these sets in the following case-by-case discussions: 

1. (zq r) G 4*1. According to Theorem 12, any (z/, r) G 4q is c-r unstable. Thus part F 
is true over 4q if one can show that (u, r) ^ r o if (u,t) G 4q. Because (u,r) £ r o 
if z^ 2 > 1 (see Eq. (4.172)), the proof for case 1 is completed if one can show that 
(z r) r o if r < v 2 < 1. 

To proceed, note that Eq. (4.175) and the facts that r o (0) = 0 and r c ( 1) = 1 imply 
that 

zA < r Q (u 2 ), v 2 < 1 (4.189) 

Thus r < T 0 {y 2 ) if r < v 2 < 1. As a result of Eq. (4.172), this in turn implies that 
(u, r) ^ r o if r < v 2 < 1. As such part F is true over 4q. 

2. (u, t ) G $ 2 . According to Theorem 16, any (zq r) G $2 is c-r unstable. Also, according 

to Eq. (4.172), (u,t) (Jl r o if (u,t) G $2. Thus part F is true over $2. 

3. (zy r) G 4*3. According to Theorem 18, any (za r) G $3 is c-r stable. Because r o (0) = 0, 
Eq. (4.172) implies that (u, r) G r o if (zq r) G 4*3. Thus part F is true over $3. 

4. (zq r) G 4*4. According to Theorem 18, any (u, r) G 4*4 is c-r stable. Because r 0 (l) = 1, 
Eq. (4.172) implies that (u, r) G r o if (u, r) G 4*4. Thus part F is true over 4*4. 

5. (zq r) G 4*5. According to Theorem 18, any (is, r) G 4*5 is c-r stable. On the other 
hand, Eqs. (4.175) implies that 

t 0 {v 2 ) < \fv 2 . 0 < v 2 < 1 (4.190) 

i.e., r G ( u 2 ) < y/v 2 = |z/| if 0 < v 2 < 1. This coupled with Eq. (4.172) implies that 
(z/, r) G r o if (z/, r) G $5. Thus part F is true over 4*5. 

6. (u, r) G 4^. For this case, we have (i) 0 < v 2 < C2, and (ii) v 2 < t < Vu 2 . To 
proceed, Note that Eqs. (4.160)-(4.162) imply that 

I+{v 2 ) < K+{ v 2 ), 0 < v 2 < c 2 (4.191) 

v 2 < L-{v 2 ) < J+{v 2 ) < VzA, 0 < v 2 < c 2 (4.192) 

and 

/ + (zA) < L-{v 2 ) < J + ( zA) < L+(z/ 2 ), 0 < z' 2 < c 2 (4.193) 

Because Eq. (4.191) contradicts Eq. (4.132), Eq. (4.132) cannot be satisfied by any 
(zy r) G 4^. Moreover, by using Eq. (4.192), it can be shown that 

g/W = 4L 1 ’ 1 ) U 4^_ 1,2 ^ U (4.194) 
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(4.195) 

(4.196) 


where \p(_ 1,2 ) . and are nonempty disjoint sets defined by 

^,( 1 . 1 ) drf |^ )T j|o < p 2 < c 2 and is 2 < r < L-(v 2 )} 

^,( 1 , 2 ) def ^ (is, 7 ~ ) 1 0 < v 2 < c 2 and L-(u 2 ) <r< J+(is 2 )} 

and 

^L 1,3) = f {(zqr) |0 < v 2 < c 2 and J+(is 2 ) < r < Vis 2 } (4.197) 

Thus any (is, r) G \P * must fall into one and only one of the following three sub-cases: 
(i) (is,r) G ^L 1,:l) , (ii) (is,r) G tL 1,2) , and (iii) (v,r) G 

Let (is,r) E By using the relation L_(z/ 2 ) < J + (is 2 ) which follows from 

Eq. (4.192) or Eq. (4.193), it is seen that Eq. (4.131) cannot be true for the current 
sub-case where v 2 < r < L-(is 2 ). Also, the second part of Eq. (4.133), i.e., L-(is 2 )< 
r < L + (u 2 ). cannot be true for the sub-case. Moreover, for a reason given earlier, 
Eq. (4.132) also cannot be true for the sub-case. According to part A of Theorem 32, 
the above results imply that any (v, r) G tP^ 1 ’ 1 is c-r unstable. On the other hand, 
because r 0 (v 2 ) = L-(v 2 ) if 0 < v 2 < (see Eqs. (4.156) and (4.171)), one concludes 
that r < r Q (u 2 ) and thus (is, r) T 0 if ( is , r) G tP^’ 1 ). As such it has been shown that 
part F is true over tP^"’ 1 ). 

Let (is,r) E tL 1,2) - It follows from Eq. (4.193) that Eq. (4.133) is satisfied by any 
(is, t) with L- ( is 2 ) < r < J + (is 2 ). According to part A of Theorem 32 and Eq. (4.196), 
this implies that any (is, r) in the current sub-case is c-r stable. On the other hand, 
because r 0 (is 2 ) = L_(is 2 ) if 0 < v 2 < C 2 , one concludes that r > r 0 (is 2 ) and thus 
(is,r) E T 0 if (is,r) E \P^’ 2 \ As such, it has been shown that part F is true over 
\p(b 2 ). 

Let (is, r) E \P^’ 3 \ Obviously Eq. (4.131) is true for the current sub-case where 
J + (v 2 ) < r < According to part A of Theorem 32, this implies that any (is, r) 

in the current sub-case is c-r stable. On the other hand, because (i) r Q (is 2 ) = L-(v 2 ) 
if 0 < v 2 < C 2 , and (ii) the relation L-(u 2 ) < J + (is 2 ) is a part of Eq. (4.193), one 

concludes that r > r 0 (is 2 ) and thus (is, r) G T 0 if (is, r) G \P^’ 3 ^. As such, it has been 

(1 3 ) 

shown that part F is true over \P_ . 

It has been shown that part F is true over each of the three nonempty disjoint sets 
tP^’ 1 ), \P^’ 2 ' ) , and \P^ L,3 - ) . Eq. (4.194) now implies that part F is true over \P^. 

7. (is, r) E \P^. For this case, we have (i) v 2 = C 2 , and (ii) is 2 < r < Vu 2 . To proceed, 
Note that Eqs. (4.163) implies that 

< I+(is 2 ) = K + (is 2 ) = L_(is 2 ) < J+(is 2 ) < Vis 2 < L+(u 2 ), is 2 = c 2 (4.198) 

With the aid of Eq. (4.198), it can be shown that 

ip^) = vp^ 2,1 ) u vp^ 2,2 -* U vp^ 2,3 ^ U tP^ 2,4 ^ (4.199) 
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(4.200) 


where T^ 2 ’ 3 ), and V I' ( J' 4) are nonempty disjoint sets defined by 

^(24) def {( u , t )\ u 2 = C2 and t/ 2 < T < L_(i/ 2 )} 

^(2,2) def t)\u 2 = C 2 and r = L_(^ 2 )} (4.201) 

^(2,3) def ^ u ^ t ^ u 2 _ C2 an( j 2) < r < J + (z/ 2 )} (4.202) 

and 

^(2,4) def r )|j,2 _ C2 and J + ^ 2 j <r< Vu 2 } (4.203) 

( 2 ) 

Thus any (is, r) G T _ must fall into one and only one of the following four sub-cases: 

(i) (u,t) G (ii) (u,t) G ^ 2,2 \ (iii) (u,t) G and (iv) (u, r) G 

Let (u,t) G By using the relation L_(z/ 2 ) < J + (u 2 ) which follows from 

Eq. (4.198), it is seen that Eq. (4.131) cannot be true for the current sub-case where 
v 2 < r < L-(u 2 ). Moreover, by using the relation I+{y 2 ) = K+{y 2 ) = L-(u 2 ) which 
also follows from Eq. (4.198), it is seen that Eq. (4.132) also cannot be true for the 
sub-case. In addition, the second part of Eq. (4.133) also cannot be true for the sub- 
case. According to part A of Theorem 32, this implies that any (u, r) G is c-r 

unstable. On the other hand, because r 0 (u 2 ) = L-(v 2 ) if v 2 = C 2 , one concludes that 
t < r 0 ( v 2 ) and thus (u, r) ^ T 0 if (u, r) G W^ 2,1 ^. As such it has been shown that part 
F is true over ^^ 2,1 \ 

Let (v, t) G ^ 2,2 \ By using the relation I + (v 2 ) = K + (u 2 ) = L_(i/ 2 ) which follows 
from Eq. (4.198), it is seen that Eq. (4.132) is true for the current sub-case where 
r = L-(u 2 ). According to part A of Theorem 32, this implies that any (v. r) G 
is c-r stable. On the other hand, because t 0 {v 2 ) = L_( v 2 ) if v 2 = C2, one concludes 

that r = t 0 (u 2 ) and thus (v, r) G T 0 if (u, t ) G T f _ 2 ' 2) • As such it has been shown that 

(2 2 ) 

part F is true over \E , _’ . 

Let (u,t) G fb^ 2 ’ 3 ). By using the relation I + (v 2 ) = L-{v 2 ) < J + (u 2 ) < L + {v 2 ) 
which follows from Eq. (4.198), it is seen that Eq. (4.133) is true for the current case 
where L_(i/ 2 ) < r < J+( u 2 ). According to part A of Theorem 32, this implies that 
any (u,t) G \b^ 2 ’ 3 ' ) is c-r stable. On the other hand, because r 0 (u 2 ) = L-(v 2 ) if 
v 2 = C2, one concludes that r > r 0 {y 2 ) and thus (u, r) G T 0 if (u, r) G W^ 2,3 ^. As such 
it has been shown that part F is true over \E_ ’ . 

Let (14 r) G \b^ 2 ’ 4 \ Obviously Eq. (4.131) is true for the current sub-case where 
J+(u 2 ) <r< Vu 2 . According to part A of Theorem 32, this implies that any (u, r) in 
the current sub-case is c-r stable. On the other hand, because (i) t 0 {v 2 ) = L _ ( v 2 ) if 
v 2 = C2, and (ii) the relation L_(v 2 )< J_i_(i/ 2 ) is a part of Eq. (4.198), one concludes 
that r > t 0 {v 2 ) and thus (u,t) G T 0 if (u,t) G 'J>^ 2,4 \ As such, it has been shown 
that part F is true over W^ 2 ’ 4 . 
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It has been shown that part F is true over each of the four nonempty disjoint sets 
ih^ 2 ’ 1 ) , ^ 2,2 \ d/^ 2 ’ 3 ), and d/^ 2,4 ). Eq. (4.199) now implies that part F is true over 

8. (u,t) G \F^. For this case, we have (i) C2 < v 2 < C3, and (ii) u 2 < r < Vv 2 . To 

proceed, Note that Eqs. (4.164) implies that 

v 2 < K+(v 2 ) < L-(u 2 ) < I+(v 2 ) < J+(v 2 ) < < L + (u 2 ), c 2 < v 2 < c 3 

(4.204) 

With the aid of Eq. (4.204), it can be shown that 

_ vi/i 3 - 1 ) y vi/C 3 . 2 ) y g>(3,3) y ^(3,4) (4.205) 

f3 1) (3 2 ) (3 3 ) (3 4 ) 

where 'LL ’ , \F_ , \F_ , and 'LL ’ ' are nonempty disjoint sets defined by 

^,(3,1) def _r^ r )| C2 < jy 2 < C3 an d v 2 <r < K + (v 2 )} (4.206) 

^,(3,2) def j^ r )| C2 < p 2 < C3 and K + (v 2 ) < r < I + (u 2 )} (4.207) 

^(3,3) def ^^ r ^| C2 < v 2 < C3 and I+{y 2 ) <r < J+{y 2 )} (4.208) 

and 

Ijy(3,4) def _r^ r )| C2 < v 2 < C3 and J + (^ 2 ) < T < Vu 2 } (4.209) 

/o\ 

Thus any (u,r) G \F_ must fall into one and only one of the following four sub-cases: 
(i) (zq r) G T^ 3 ’ 1 ^, (ii) (zq r) G T^ 3 ’ 2 ^, (iii) (zq r) G T^ 3 ’ 3 ^, and (iv) (zq r) G T^ 3,4 ^. 

Let (z/, t) G By using the relation Ay_(zA) < J + (z/ 2 ) which follows from 

Eq. (4.204), it is seen that Eq. (4.131) cannot be true for the current sub-case where 
v 2 <r < K + (v 2 ). Moreover, obviously Eq. (4.132) is also not true for the sub-case. 
In addition, by using the relation K + (i/ 2 ) < L_(z/ 2 ) < I + {u 2 ) which also follows from 
Eq. (4.204), one concludes that Eq. (4.133) also can not be true for the sub-case. 
According to part A of Theorem 32, the above results imply that any (zq r) G T^ 3,1 ^ 
is c-r unstable. On the other hand, because t 0 (u 2 ) = K+(v 2 ) if c 2 < i/ 2 < c 3 , one 
concludes that r < t 0 {v 2 ) and thus (zq r) ^ T 0 if (zq r) G T^ 3,1 ^. As such it has been 
shown that part F is true over T_’ . 

(3 2 ) 

Let (u,t) G \F_ . Obviously Eq. (4.132) is true for the current sub-case where 

K + {u 2 ) < r < I + (u 2 ). According to part A of Theorem 32, this implies that any 
(u, t ) G T^ 3,2 '* is c-r stable. On the other hand, because r 0 (u 2 ) = K + (u 2 ) if c 2 < v 2 < 
c 3 , one concludes that r > r 0 ( z/ 2 ) and thus (zq r) G r o if {u. r) G T^ 3,2 ^. As such it 
has been shown that part F is true over T_ ’ . 

Let (zq t) G T^ 3,3 ^. By using the relation L_( u 2 ) < /_|_(z/ 2 ) < J+(z^ 2 ) < L + (u 2 ) 
which follows from Eq. (4.204), it is seen that Eq. (4.133) is true for the current case 


NASA/TM— 2005-2 1 3627 


35 



where I+(v 2 ) < r < J+(v 2 ). According to part A of Theorem 32, this implies that 
any (u,t) G is c-r stable. On the other hand, because (i) t 0 {v 2 ) = K + (u 2 ) 

if C2 < v 2 < C3, and (ii) the relation K + (i' 2 ) < I+{y 2 ) is a part of Eq. (4.204), one 

concludes that r > r 0 (u 2 ) and thus (u,t) G T 0 if (u,t) G As such it has been 

(3 3) 

shown that part F is true over \E'_’ . 

Let (u,t) G Obviously Eq. (4.131) is true for the current sub-case where 

J+(u 2 ) < r < \fu 2 . According to part A of Theorem 32, this implies that any (u, r) 

in the current sub-case is c-r stable. On the other hand, because (i) T 0 {y 2 ) = K + {y 2 ) 
if C2 < v 2 < C3, and (ii) the relation K+(v 2 ) < J + (u 2 ) is a part of Eq. (4.204), one 
concludes that r > t 0 (v 2 ) and thus (u, r) G r o if (u, r) G As such, it has been 

shown that part F is true over 

It has been shown that part F is true over each of the four nonempty disjoint sets 
v]/! 3 ’ 3 )^ a nd v^O’ 4 ) gq (4.205) now implies that part F is true over 4/^. 

9. (iq r) G For this case, we have (i) C3 < v 2 < 1, and (ii) v 2 < r < \f\P- . To 

proceed, Note that Eqs. (4.165)-(4.168) implies that 

v 2 <K + (v 2 )<V^<J + (p 2 )<I+(v 2 )<L + (v 2 ), c 3 <u 2 < 1 (4.210) 

With the aid of Eq. (4.210), it can be shown that 

_ ^(44) y ^(4,2) (4.211) 

where and T 1 ) 4 ' 2 ' 1 are nonempty disjoint sets dehned by 

^(_ 4 d) def {(V ?r )| C3 < < 1 and v 2 <T< K+(u 2 )} (4.212) 

and 

^(4,2) def < v 2 < 1 an d K + (v 2 ) < r < Vu 2 } (4.213) 

Thus any (u, r) G must fall into one and only one of the following two sub-cases: 
(i) (u, t) G T^ 4,1 ' 1 and (ii) (1/, r) G \t'^ 4,2 \ 

Let (u,t) G ^'^ 4 ’ 1 ^. By using the relation K + (v 2 ) < J + (u 2 ) < /_|_ (j/2) which follows 
from Eq. (4.210), it is seen that none of Eqs. (4.131)-(4.133) is true for the current 
sub-case where v 2 < r < K + {v 2 ). According to part A of Theorem 32, this implies 
that any (v,t) G ’L^ 4 ’ 1 -’ is c-r unstable. On the other hand, because r 0 (u 2 ) = K + (u 2 ) 
if C3 < v 2 < 1, one concludes that r < t 0 (v 2 ) and thus (zqr) ^ T 0 if (u, r) G 
As such it has been shown that part F is true over 'L^ 4,1 \ 

Let (u, r) G ^ 4 ’ 2 . By using the relation K+(v 2 ) < \fv 2 < /_ |_(zA) which follows 
from Eq. (4.210), it is seen that Eq. (4.132) is true for the current sub-case where 
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K + (u 2 ) < r < y/JS. According to part A of Theorem 32, this implies that any 
(u, t) G is c-t stable. On the other hand, because t 0 (v 2 ) = K + (u 2 ) if C3 < v 2 < 

1, one concludes that r > r G ( v 2 ) and thus (u, r) G T 0 if (u, r) G As such it has 

been shown that part F is true over F*2 2 \ 

It has been shown that part F is true over each of the two nonempty disjoint sets 
W^ 4,1) and Eq. (4.211) now implies that part F is true over 

10. (u, t) G F+f For this case, we have (i) 0 < v 2 < C 3 , and (ii) r > \fv 2 . To proceed, 
Note that Eqs. (4.160)-(4.165) imply that 

I + (is 2 ) < \4a, 0 < v 2 < c 3 (4.214) 

K + {v 2 ) < Vu 2 , 0 < v 2 < c 3 (4.215) 

and 

L_(v 2 ) < Vu 2 , 0 < v 2 < C 3 (4.216) 

By using Eqs. (4.214) and (4.215), one concludes that Eq. (4.135) is true for the current 
case where r > \f~J 2 . According to part B of Theorem 32, this implies that any (zy r) 
in the current case is c-t stable. On the other hand, because (i) t 0 (v 2 ) = L-(v 2 ) if 
0 < v 2 < C 2 , and (ii) r 0 {v 2 ) = K+(v 2 ) if C 2 < v 2 < C 3 , Eqs. (4.215) and (4.216) imply 
that t > T 0 {y 2 ) and thus (v, r) G r o if (u,r) G F+f As such, it has been shown that 
part F is true over F+\ 

11. (u,t) G F^. For this case, we have (i) C 3 < v 2 < 1, and (ii) r > \fv 2 . To proceed, 
Note that Eqs. (4.166)-(4.168) imply that 

K+(u 2 ) < Vu 2 < J+(u 2 ) < I+(u 2 ), c 3 < v 2 < 1 (4.217) 

and 

L_(i/ 2 ) < J + (i/ 2 ) < i+(i/ 2 ) < T_|_(z^ 2 ), C 3 < v 2 < 1 (4.218) 

By using Eq. (4.217), one has 

^( 2 ) = ^( 2 , 1 ) u ^( 2 , 2 ) u \J/(2>3) (4.219) 

(2 1 ) (2 2 ) (2 3 ) 

where F£’ x ', F+’, and F^’ J; are nonempty disjoint sets defined by 

F^’i) = f {(u,t)\cs < v 2 < 1 and Vv 2 < r < J + (p 2 )} (4.220) 

^( 2 , 2 ) def |^^ r ^| C3 < u 2 < 1 and J + {v 2 ) <r < I + {v 2 )} (4.221) 

and 

^(2,3) def |^ r )| C3 < ^ < i and r > / + (z/ 2 )} (4.222) 
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( 2 ) 

Thus any (u,t) E must fall into one and only one of the following three sub-cases: 
(i) (v,t) E (ii) (u, r) E 'L+’ 2) , and (iii) (u,r) E 'L+’ 3) . 

Let (v,t) E ’ . Eq. (4.134) is true for any (v,t) in the current sub-case where 
Vu 2 < t < J- \-{y 2 ). According to part B of Theorem 32, this implies that the any 
(u, r) E is c-r stable. On the other hand, because (i) t 0 {v 2 ) = K + (p 2 ) if 

C 3 < r < 1, and (ii) the relation K + (u 2 ) < Vv 2 is a part of Eq. (4.217), one concludes 
that r > r 0 (u 2 ) and thus (u, r) E T 0 if (u, r) E As such it has been shown that 

part F is true over . 

Let (zq r) E T^ 2,2 ^. By using Eq. (4.218), one concludes that Eq. (4.136) is true 
for the current case where J+(v 2 ) < r < f+(V 2 ). According to part B of Theorem 
32, this implies that any (zq r) E ' is c-r stable. On the other hand, because (i) 
r 0 (u 2 ) = K + (is 2 ) if C 3 < v 2 < 1, and (ii) the relation K + (v 2 ) < J + ( v 2 ) is a part of 

Eq. (4.217), one concludes that r > r 0 (u 2 ) and thus (zq r) E T 0 if (zq r) E As 

(2 2 ) 

such it has been shown that part F is true over . 

Let (z/, r) E T^’ 3) . By using the relation K + (u 2 ) < J_|_(zA) which follows from 
Eq. (4.217), one concludes that Eq. (4.135) is true for the current sub-case where 
r > I + (u 2 ). According to part B of Theorem 32, this implies that any (u. r) E Tf, 2 ' 3 * 
is c-r stable. On the other hand, because (i) t 0 (u 2 ) = K + (u 2 ) if C 3 < v 2 < 1 , and (ii) 
the relation K + {v 2 ) < /+(z^ 2 ) is a part of Eq. (4.217), one concludes that r > t 0 (u 2 ) 
and thus (u. r) E T 0 if (u, r) E 41+’ . As such, it has been shown that part F is true 

VTr( 2 > 3 ) 

over W^|_ . 

It has been shown that part F is true over each of the three nonempty disjoint sets 
^(24), ^ an( j gjq. (4.219) now implies that part F is true over 41^. 

It has been established that part F is true over each of the sets mentioned in the 
paragraph immediately following Eq. (4.188). Because any (z/, r) must belong to one and 
only one of these sets, the proof of part F is completed. 

Finally, with the aid of Theorems 4 and 16, one can obtain part G from part F. QED. 

As promised earlier, a proof for Theorem 34 will be provided in the remainder of the 
paper. As a preliminary, we have Theorem 36. 

Theorem 36. In the domain 0 < x < 1, (A) I+(x), J+(x), K + (x), and L_(x) are 
strictly monotonically increasing while L+(x) is strictly monotonically decreasing. More- 
over, we have (B) 


co 

V 

£ 

)>o, 

0 < x < 1 

(4.223) 

3 > J + (x 

)>0, 

0 < x < 1 

(4.224) 

1 > K + (x ) > 0, 

0 < : 

c < 1 

(4.225) 

3 > L_(x 

) >0, 

0 < a 

; < 1 

(4.226) 
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L+(x) > 3, 


0 < x < 1 


(4.227) 


Proof. Let f'(x) = df(x)/dx for any function / of x. Then (i) Eqs. (4.91) and (4.98) 
imply that 


j , v 3x — 2 + 2\/3x 2 — 3x + 1 

4* = , „ = — > o, 


'3x 2 — 3x + 1 


0 < x < 1 


(4.228) 


(ii) Eqs. (4.92) and (4.100) imply that 


= 


—x s + 6x 2 — x + 2 + 4y / 2(x 3 — x + 2) 

(2 — x) 2 y2(x^^aT+~2) 

x 2 (l — x) + 5a: 2 + 1 + (1 — x) + 4i,/2(x 3 — a: + 2) 
(2 — x) 2 y/2(x 3 — x + 2) 


(iii) Eqs. (4.114) and (4.118) imply that 

= Vl-2x + 5x 2 -(l-x) 

+ 2a; 2 \/l — 2x + 5a: 2 

(iv) Eqs. (4.122) and (4.126) imply that 

2 U — x — 2^2(2 — a- — a; 2 ) 

(a;) = — ^ J- > 0, 

x 2 -\/2(2 — a; — a; 2 ) 

and (v) Eqs. (4.122) and (4.126) imply that 

2 [4 - a: + 2^2(2 - x - x 2 ) 

L'Ax) = ^ < 0, 

x 2 a/2(2 — a; — a; 2 ) 


(4.229) 


> 0, 0 < x < 1 


0 < x < 1 


(4.230) 


0 < x < 1 


(4.231) 


0 < x < 1 


(4.232) 


Thus part A is true. 

Moreover, by using (i) Eqs. (4.91), (4.92), (4.114), and (4.122), and (ii) L’hopital’s 
rule, one has (i) 


lim i+(x) = lim J+(x) = lim L_(x) = lim L + (x) =3, and lim K + (x) = 1 

X^l~ X^l~ X^l~ x^l~ x^l~ 


(4.233) 


lim I+(x) = lim ( 3 -\ — ®x 3 \ _ (—3) = 0 

c -o+ +v ' x^o+ V V3x 2 - 3x + 1 J V ' 
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lim J+(x) = 0 

x^0+ 


(4.235) 


and (v) 


1 / 5r — 1 \ 1 

lim K+(x) = lim - 1 + , = _(i - 1) = 0 

a;— >-o+ x >o~i" 2 \ Vl — 2x + 5a; 2 J 2 


lim L-(x ) = lim —1 H 2(1 + 2a:) =—1 + 1 = 0 

a;^o+ \/-{‘2 — ./• .r 2 ) 


(4.236) 


(4.237) 


part B now follows from Part A and Eqs. (4.233)-(4.237). QED 

An immediate result of Theorem 36 and the fact that 0 < x < y/x < 1 if 0 < x < 1 is 
given in Theorem 37. 

Theorem 37. We have 


x < \[x < L + (x), I+(x ) < L + (x), J+(x ) < L + (x), 
K + (x) < L + (x ) and L-(x) < L + (x), 0 < x < 1 


(4.238) 


Theorem 37 is but one of many algebraic relations that are needed in the proof of 
Theorem 34. Note that, in establishing other needed relations, several inequalities that 
involve the four prinicipal square roots that appear in the definitions of I±(x), J±(x), 
K±(x), and L±(x), i.e., 


'3a: 2 — 3a: + 1 > 0, 


— oo < x < +oo 


(4.239) 


^2(x 3 - x + 2) > 0, 
'l - 2x + 5a; 2 > 0, 


0 < x < 2 


— oo < x < +oo 


(4.240) 

(4.241) 


/ 2(2-a:-a: 2 ) > 0, 


-2 < x < 1 


(4.242) 


(which follow from Eqs. (4.97), (4.100), (4.117), and (4.124), respectively) will be used 
repeatedly. Also to be used often is the following algebraic property: 

Property I. Let a > 0 and b > 0. Then 

f>0 a — b > 0 


a 2 —b 2 { = 0 o-6 = 0 


<0 a — b < 0 


(4.243) 
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With the above preparations, a set of relations will be given in Theorems 38-48. 

Theorem 38. We have 


( > 0 if 0 < x < 3 — 2^/2 


x — I+(x) 


0 if x = 3- 2V2 


l < 0 if 3 - 2V2 <x<l 


(4.244) 


Proof. Let 0 < x < 1 throughout the proof. Then Eq. (4.91) implies that 


x — I+(x) = 


X 


2 - 3a: + 2 - 2V3x 2 - 3a: + 1 


x 


(4.245) 


With the aid of Property I, Eq. (4.244) is a result of Eq. (4.245) and the following relations: 
(i) Eq. (4.239); (ii) 

x 2 -3x + 2 = (x- l)(x-2) > 0 (4.246) 


(hi) 

(x 2 — 3x + 2) 2 — ^2\/3x 2 — 3x + l^j = x 2 (x 2 — 6x + 1) 

= x 2 \x - (3 + 2V2)] \x - (3 - 2V2) 

and (iv) 0 < 3 - 2\[2 < 1 < 3 + 2y/2. QED. 

Theorem 39. We have 


(4.247) 


x < K + (x), 


0 < x < 1 


(4.248) 


Proof. Let 0 < x < 1 throughout the proof. Then Eq. (4.114) implies that 

y/l — 2x + 5x 2 — (2x 2 — x + 1) 


K- |_(x) — x = 


2x 


(4.249) 


With the aid of Property I, Eq. (4.248) is a result of Eq. (4.249) and the following relations: 
(i) Eq. (4.241); (ii) 

2x 2 — x + 1 = 2(x — 1/4) 2 + 7/8 > 7/8 (4.250) 


and (iii) 

QED. 


1 — 2x + 5x 2 ) — (2x 2 — x + l) 2 = 4x 3 (1 — x) > 0 


(4.251) 
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Theorem 40. Let C 3 be the constant defined in Eq. (4.157). Then 

I > 0 if 0 < x < c 3 

= 0 if x = C 3 
<0 if c 3 < x < 1 

Proof. Let 0 < x < 1 throughout the proof. Then Eq. (4.91) implies that 

x\fx — 3x + 2 — 2\/3x 2 — 3x + 1 


(4.252) 


x — I+(x) = 


x 


(4.253) 


With the aid of Property I, Eq. (4.252) is a result of Eq. (4.253) and the following relations: 
(i) Eq. (4.239); (ii) 


x 


\fx — 3x T 2 — ^1 — \Zx) [l T 2y/x T (1 — x) !> 0 


(4.254) 


(in) 


(x\[x — 3x + 2 ) 2 — (2\/3x 2 — 3x + l ) 2 = x 3 — 6x 5 ^ 2 — 3x 2 + 4a : 3 / 2 


= x 3 ! 2 (\fx + l) ^\/x — 


7 + V33 


x — 


7-V33\ 

2 J 


(4.255) 


(iv) 0< (7- V33)/2< 1 < (7 + V33)/2; and (v) c 3 = (7 - ^33) /2 . QED. 
Theorem 41. Let C 3 be the constant defined in Eq. (4.157). Then 

I > 0 if 0 < x < c 3 

= 0 if x = C 3 
<0 if C 3 < x < 1 

Proof. Let 0 < x < 1 throughout the proof. Then Eq. (4.92) implies that 

—Xy/x — 3x + 2^/x + 2 — ^/2(x 3 — x + 2) 


(4.256) 


\[x — j_|_ (x) = 


2 — x 


(4.257) 


With the aid of Property I, Eq. (4.256) is a result of Eq. (4.257) and the following relations: 
(i) Eq. (4.240); (ii) 


—x\fx — 3x + 2 \fx + 2 = (l — VS) (x + 3 l\[x + 2) > 0 


(4.258) 
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(iii) 


■i 2 


(-xy^ — 3x + 2 yjx + 2) — \/2[x 3 — x + 2) 
= —x 3 + 6 a : 5 / 2 + 5a : 2 — 16a: 3 / 2 — 6 a: + 8 \fx 


n , o M r ^\l r 7 + V^\ { ^ 7 - VS3\ 

= y/x{2 — x) (v x + 1 j ( v ■ x - J ( V ■ x 2 j 


(4.259) 


(iv) 0 < (7 - V33) /2 < 1 < (7 + y/33)/2; and (v) c 3 = (7 - y/33)/2 

Theorem 42. We have 


. QED. 


K+(x) < yfc, 


0 < x < 1 


(4.260) 


Proof. Let 0 < x < 1 throughout the proof. Then Eq. (4.114) implies that 


\[x — K + (x ) = 


2xy/x — x + 1 — yl — 2 a: + 5 a : 2 


2 a: 


(4.261) 


With the aid of Property I, Eq. (4.260) is a result of Eq. (4.261) and the following relations: 
(i) Eq. (4.241); (ii) 

2xyfx — x + 1 = 2xy/x + (1 — x) > 0 (4.262) 


and (iii) 


[2x\[x — x + l) — 1 — 2 a: + 5a : 2 J = 4x\Jx(l — x) (l — \fx) > 0 

QED. 

Theorem 43. Let C 4 be the constant defined in Eq. (4.158). Then 

' > 0 if 0 < x < C 4 
fx — L_(x) ^ =0 if x = C 4 

<0 if C 4 < x < 1 


(4.263) 


(4.264) 


Proof. Unless specified otherwise. Let 0 < x < 1 in this proof. Then Eq. (4.122) 
implies that 


x — L-(x) 


2^/2(2 — x — a: 2 ) — (4 — x — xy/x) 


x 


(4.265) 
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To proceed, note that 



2 ^ 2(2 - x - x 2 ) 

2 2 
— (4 — x — x\fx) — — x\Jx g{x) 

(4.266) 

where 

5 f(x) c = + 2 x + 9y^ — 8 , x > 0 

(4.267) 

Because (i) 




g'(x) = Sy/x/2 + 2 + 9/(2y/x) = 3/(2^) [{Vx + 2/3 ) 2 + 23/9] > 0, 

x > 0 (4.268) 

and (ii) 

P( 0 ) = 

= —8 and g(l) = 4 

(4.269) 


one concludes that g(x ) is strictly monotonically increasing in the interval 0 < x < 1 and 
there is one and only one real root of g(x) = 0 in this interval. By using the standard 
formula for the roots of a cubic equation, it can be shown that this root is given by x = C 4 . 
Moreover, Eqs. (4.268) and (4.269) imply that: (i) g(x ) < 0 if 0 < x < c\\ (ii) g(x) = 0 if 
x = C4; and (iii) g(x) > 0 if C4 < x < 1. As such Eq. (4.266) implies that 

>0 if 0 < x < C 4 


= 0 if x = C 4 (4.270) 

<0 if C 4 < x < 1 

With the aid of Property I, Eq. (4.264) is a result of Eqs. (4.265) and (4.270), and the the 
following relations: (i) Eq. (4.242); and (ii) 

4 — x — x\[x = 2 + (1 — x) + (l — x\fx ) > 2 , 0 < x < 1 (4.271) 


2^/2(2 — x — x 2 ) — (4 — x — x\fx)‘ 


QED. 

Theorem 44. We have 


L_(x) < J + (x ), 


0 < x < 1 


(4.272) 


Proof. Let 0 < x < 1 throughout this proof. Then Eqs. (4.92) and (4.122) imply that 


J + (x) - L-{x) = - ^ + 2) + 2(2 - *K/2(2 - x - - (8 - lx - 2^) (4273) 

Xy^i X ) 


Let 


P(x) = f x^/2(x 3 — x + 2) + 2(2 — x)y/2{2 — x — x 2 ) + (8 — 4x — 2x 2 ) (4.274) 
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and 


(4.275) 


Then (i) 


/ 3±(x ) *== 4 \J (a; 3 — x + 2)(2 + x) ± \/l — x (8 + 3a: — x 2 ) 


[«/+(») ~L-(x)]P(x) 

8x(2 — x)^/(x3 — x + 2)(1 — a:) (2 + x) + 2a; 5 — 12a; 4 + 6a; 3 + 36a; 2 — 32a: 

x(2 — x) 

8x{2 — x)a/(3; 3 — x + 2)(1 — x)(2 + x) — 2x{l — x){2 — x)(8 + 3x — x 2 ) 

x(2 — x) 

= 2\/l — x /3- (x) 
and (ii) 

/ 3 _ (a;)/?_|_(x) = a; 5 + 9a; 4 + 31a; 3 + 39x 2 + 16a; 

Thus 


[.J + (x) — L-(x)\(3(x)(3+(x) = 2^1 — x(x 5 + 9a; 4 + 31a; 3 + 39a; 2 + 16a;) 
By using (i) Eqs. (4.240) and (4.242); (ii) \J (a; 3 — x + 2)(2 + x) > 0; (iii) 

8 - 4a: - 2a; 2 = 2 + 4(1 - x) + 2(1 - x 2 ) > 2 


and (iv) 

8 + 3a; — x 2 = 7 + 3a; + (1 — x 2 ) > 7 
it follows from Eqs. (4.274) and (4.275) that 

/3(x) > 0 and P+(x) >0, 0 < x < 1 

Eq. (4.272) is a result of Eq. (4.281) and the fact that the expression on the right 
Eq. (4.278) is positive everywhere in the interval 0 < x < 1. QED 

Theorem 45. We have 


( > 0 if 0 < a: < 3/11 


K+(x)-I+(x) < =0 if x = 3/11 


l < 0 if 3/11 < x < 1 


(4.276) 

(4.277) 

(4.278) 

(4.279) 

(4.280) 

(4.281) 
side of 

(4.282) 
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Proof. Eqs. (4.91) and (4.114) imply that 


, , (3 - 5a;) - (4\/3a; 2 - 3a; + 1 - Vl - 2a; + 5a; 2 ) 

K + (x) — I+(x) = 0 < x < 1 (4.283) 

2.r 

By using (i) Eqs. (4.239) and (4.241), and (ii) 

(±y/Zx 2 - 3a;+ l) 2 - (Vl - 2a; + 5a; 2 ) 2 = 43a; 2 - 46a; + 15 

= 43 [(a; - 23/43) 2 + 116/(43) 2 ] > 116/43, -oo < x < oo 
an application of Property I leads to the conclusion 

4 \J 3a; 2 — 3a; + 1 — \[ 1 — 2a; + 5a; 2 > 0, — oo < x < +oo 

Moreover, we have 


<0 if a; > 3/5 


3 -5a; 


>0 if x < 3/5 
Combining Eqs. (4.283), (4.285) and (4.286), one has 

K + (x) — I+(x) < 0, 3/5 < x < 1 


(4.284) 


(4.285) 


(4.286) 


(4.287) 


To proceed, let 

£(a:) d = ^3 - 5a; + 4\/3a; 2 - 3a; + 1 - \f\ - 2a; + 5a; 2 ) , 0 < x < 1 (4.288) 

and 

£±(x) d = 2y/ (3x 2 - 3a; + 1)(1 - 2a; + 5a; 2 ) ± (7a; 2 - 5x + 2), 0 < x < 1 (4.289) 

Then Eqs. (4.285) and (4.286) imply that 

£(a;) >0, 0 < x < 3/5 (4.290) 

In addition, by using (i) \/(3a' 2 — 3a; + 1)(1 — 2a' + 5a; 2 ) > 0, — oo<a;<+oo (which follows 
from Eqs. (4.239) and (4.241)); and (ii) 

7a; 2 -5a; + 2 = 7 [(a; - 5/14) 2 + 31/196] >31/28, -oo < x < +oo (4.291) 
one has 

£+(x) >0, 0 < x < 1 (4.292) 
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Combining Eqs. (4.290) and (4.292), one arrives at the conclusion: 


£(x) £+(x) >0, 0 < x < 3/5 

Next, Eqs. (4.283), (4.288), and (4.289) imply that (i) 

[K + (x) — I + (x)]£(x) = — — , 0 < x < 1 

x 

and (ii) 

£-(x) £+(x) = llx 4 — 14a: 3 + 3a: 2 = lla: 2 (a: — l)(ar — 3/11), 0 < x < 1 

Thus 


[K + (x) — I+(x)] £(x)£+(x) = llx(x — l)(a: — 3/11), 0 < x < 1 

It follows from Eqs. (4.293) and (4.296) that 

' > 0 if 0 < x < 3/11 
K + (x) — I + (x) < =0 if x = 3/11 

, < 0 if 3/11 < x < 3/5 

Eq. (4.282) is an immediate result of Eqs. (4.287) and (4.297). QED. 
Theorem 46. We have 

(>0 if 0 < a: < 3/11 
L_(x) — I+(x) < =0 if x = 3/11 

, < 0 if 3/11 < x < 1 


Proof. Eqs. (4.91) and (4.122) imply that 
L-(x) — I + (x) = — y3 — 2x — y / 2(2 — x — a: 2 ) — \J 3a: 2 — 3a: + 1 J , 0 < x < 

By using Eq. (4.299) and the definitions 

/r(x) = f - ^3 — 2a: + \/2(2 — x — a: 2 ) + \j 3a: 2 — 3a: + 1^ , 0 < x < 1 


(4.293) 


(4.294) 


(4.295) 


(4.296) 


(4.297) 


(4.298) 


1 (4.299) 


(4.300) 


NASA/TM— 2005-2 1 3627 


47 



and 

H±{x) = f 3a; 2 — lx + 4 ± 2y/2(2 — x — a; 2 )(3a; 2 — 3a; + 1), 0 < x < 1 

one has 

[L-(x) — I+(x )] //(a;) = - — , 0 < a; < 1 

and 

H-{x) n+(x) = 33a; 4 — 42a; 3 + 9a; 2 = 33x 2 (x — l)(x — 3/11), 0 < x < 1 

In turn, Eqs. (4.302) and (4.303) imply that 

\L-(x) — I+(x)] n{x) l^+(x) = 33x(x — l)(a; — 3/11), 0 < x < 1 

By using (i) Eqs. (4.239) and (4.242); (ii) 3 — 2a;>0ifa;<3/2; and (iii) 

3a; 2 — 7x + 4 = 3(a; — l)(a; — 4/3) >0, x < 1 or x > 4/3 

Eqs. (4.300) and (4.301) imply that 

/i(x) > 0 and fj, + (x) > 0, 0 < x < 1 

Eq. (4.298) is an immediate result of Eqs. (4.304) and (4.306). QED. 
Theorem 47. We have 

f>0 if 0 < x < 3/11 or 3/11 < x < 1 
L-(x)-K+(x)\ 

1=0 if x = 3/11 


Proof. Eqs. (4.114) and (4.122) imply that 

9 — 3a; — 4^/2(2 — x — x 2 ) — \/l — 2x + 5a; 2 
L-{x) - K+(x ) = — 

By using Eq. (4.308) and the definitions 

ip(x) '= ^ ^9 — 3x + 4y/2(2 — x — a; 2 ) + \/l — 2x + 5a; 2 ^ , 

and 


0 < x < 1 


0 < x < 1 


ip ±{x ) *= f 9x 2 — 5x + 4 ± 2y/2(2 — x — a; 2 )(l — 2a; + 5x 2 ), 0 < x < 1 


(4.301) 

(4.302) 

(4.303) 

(4.304) 

(4.305) 

(4.306) 

(4.307) 

(4.308) 

(4.309) 

(4.310) 
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one has 

[L_(x) — K + (x)\ il>(x) = — ^ , 0 < a; < 1 (4.311) 

X 

and 

ip- (a;) ^(a;) = 121a; 4 — 66a; 3 + 9a; 2 = 121x 2 (x — 3/11) 2 , 0 < a; < 1 (4.312) 

In turn, Eqs. (4.311) and (4.312) imply that 

[L-(x) — K + (x)\ ip (x) ip + (x) = 121x(x — 3/ ll) 2 , 0 < x < 1 (4.313) 

By using (i) Eqs. (4.241) and (4.242); (ii) 9 — 3a; > 0 if x < 3; and (iii) 

9a; 2 -5x + 4 = 9 [(a; - 5/18) 2 + 119/324] > 119/36, -oo < x < +oo (4.314) 

Eqs. (4.309) and (4.310) imply that 

ip(x) > 0 and ip+(x) >0, 0 < x < 1 (4.315) 

Eq. (4.307) is an immediate result of Eqs. (4.313) and (4.315). QED. 


Theorem 48. Let C3 be the constant defined in Eq. (4.157). Then we have 



' > 0 

if 

0 < x < c 3 


J+(x)-I+(x) < 

= 0 

if 

x = c 3 

(4.316) 


, < 0 

if 

c 3 < x < 1 



Proof. Eqs. (4.91) and (4.92) imply that 


J+(x)-I+(x ) 


6a; 2 — 10a; + 4 


2(2 — a;)V3a; 2 — 3a; + 1 — a;>/2(a; 3 — x + 2) 
ar(2 — ar) 


0 < x < 1 

To proceed, note that Eq. (4.239) implies that 

2(2 — a;)\/3a; 2 — 3x + 1 > 0, x<2 

Also Eq. (4.240) implies that 

a;>/2(x 3 — x + 2) >0, 0 < x < 2 


(4.317) 

(4.318) 


(4.319) 
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Moreover, we have 


2(2 — a;)\/3a; 2 — 3a; + 1 


a;\/2(a; 3 — x + 2) 


= —2a; 5 + 12a; 4 — 58a; 3 + 96a; 2 — 64a; + 16 
= 2(1 — x)(x 4 — 5a; 3 + 24a; 2 — 24a; + 8) 

= 2(1 — x) \x 2 {x 2 — 5a; + 6) + 2(9a; 2 — 12a; + 4)] 

= 2(1 - x) [x 2 (x -2)(x-3) + 2 (3a; - 2) 2 ] >0, 0 < x < 1 


(4.320) 


With the aid of Eqs. (4.318)-(4.320), an application of Property I leads to the conclusion 


2(2 - x) ^3a; 2 - 3ar + 1 - xy/2[rf - x + 1) > 0, 0 < x < 1 (4.321) 


Next note that 

a;(2 — a;) > 0, 0 < a; < 2 


(4.322) 


6a; 2 — 10a; + 4 = 6(a; — l)(a; — 2/3) 


and 

<0 if 2/3 < x < 1 

>0 ifa;<2/3ora;>l 
By combining Eq. (4.317) with Eqs. (4.321)-(4.323), one concludes that 

J+(x) — I+(x) <0, 2/3<a;<l 

To study the case where 0 < x < 2/3, let 

rj(x) = f - 6a; 2 — 10a; + 4 + 2(2 — x) \J 3a; 2 — 3a; + 1 — a;^/2(a; 3 — x + 2) 
2 


(4.323) 


(4.324) 


0 < x < 1 
(4.325) 


and 

r}±(x) d = 2\f2^x 2 - 3a; + l)(a; 3 - x + 2) ± (-a; 3 + 10a; 2 - 9a; + 4), 0 < x < 1 (4.326) 
By using Eqs. (4.321) and (4.323), Eq. (4.325) implies that 

rj( x) >0, 0 < x < 2/3 (4.327) 


Moreover, because (i) (3a; 2 — 3a; + l)(a; 3 — x + 2) > 0, 0 < x < 2 (see Eqs. (4.239) and 

(4.240)); and (ii) 

—x 3 + 10a; 2 — 9a; + 4 = a; 2 (l — x) + (3a; — 3/2) 2 + 7/4 > 7/4, 0<x<l (4.328) 
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Eq. (4.326) implies that 

f] + (x) > 0, 0 < x < 1 (4.329) 

Next, by using Eqs. (4.157), (4.317), (4.325) and (4.326), it can be shown that (i) 

[J+(x) ~ I+(x)]i l(x) 

2\[2 x(2 — x) \J (3a; 2 — 3x + l)(x 3 — x + 2) — (a; 5 — 12a; 4 + 29a; 3 — 22x 2 + 8x ) 

x(2 — x) 

2\f2x(2 — x)\J (3a; 2 — 3a: + l)(a; 3 — x + 2) — x(2 — x)(—x 3 + ICtr 2 — 9x + 4) 

a:(2 — x) 


= V-( x ), 


0 < x < 1 


and (ii) 


(4.330) 


?7_(x)?7_|_(x) = —x 16 + 44x b — 142a; 4 + 172x 3 — 89a; 2 + 16a; = x(l — x) 6 (x z — 41x + 16) 


= a;(l — x) 3 (x — C3) ( x — 


41 + 7^33 


In turn, Eqs. (4.330) and (4.331) imply that 


0 < x < 1 


(4.331) 


[J+(x) — I+(x)] r/(x ) r/ + (x) = a;(l — x) 3 (x — c 3 ) (x — — — j , 0 < x < 1 (4.332) 


With the aid of Eqs. (4.327), (4.329) and (4.332), and the relation 

0 < c 3 < 2/3 < 1 < 41 + l' m 


(4.333) 


one concludes that 


( > 0 if 0 < x < c 3 

= 0 if x = c 3 


J+(x) — I+(x) < = 0 if x = c 3 (4.334) 

l < 0 if c 3 < x < 2/3 

Eq. (4.316) now is an immediate result of Eqs. (4.324) and (4.334). QED. 

With the above preparations, Theorem 34 can now be proved. Part A is identical to 
part A of Theorem 36. Part B can be shown using Theorems 38-48 and the two relations 


\/x < L + (x ) and I+(x ) < L + (a;), 


(4.335) 


NASA/TM— 2005-2 1 3627 


51 



which form a part of Theorem 37. Part C follows from Eqs. (4.230), (4.231), and (4.156). 
Part D was shown in Eqs. (4.233) and (4.237). QED. 

Finally, note that none of the relations 

x < J+(x), x < L_(x), and K + (x) < J+(x), 0<x<l (4.336) 

appears in Theorems 37-48. However, they can be shown using Theorem 34. As such, 
they can be considered as results of Theorems 38-48 and the relations Eq. (4.335). 
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5. Conclusions and Discussions 


With the aid of many unexpected mathematical simplifications that occur along the 
way, it has been shown in Sec. 4 that there is an explicit analytical solution to the implicit 
stability conditions stated in Theorem 3. The first and perhaps the most important “break” 
encountered is the simple relation Eq. (4.23), i.e., H(u,r,s), a quartic polynomial in s, is 
equal to the product of 4(1 — u 2 )s 2 and G ( u. r, s), a quadratic polynomial in s. Without 
Eq. (4.23) and the fortunate fact that both D(u. r, s) and F(v. r, s) are also quadratic 
polynomials in s, the relatively straightforward study of the necessary stability conditions 
Eqs. (4.25)-(4.27) (Theorem 6) as presented in Sec. 4 would have become much more 
complicated. 

Moreover, the fact that F(v, r, 1) and H(u , r, 1) can be cast into the simple factorized 
forms Eqs. (4.35) and (4.37), respectively, are instrumental in the successful effort to 
establish Eq. (4.41) as necessary conditions for stability (Theorem 12). 

With the aid of Theorems 13-15, it was shown that the special case in which (u, r) 
satisfies Eq. (4.2) and yet is c-r unstable occurs if and only if r = v 2 = 1 (Theorem 16). 
Using Theorem 16, Theorem 17 was then established to provide a set of necessary and 
sufficient stability conditions much more explicit and easier to handle than those given 
originally in Theorem 3. Based on Theorem 17, it was then shown that the c-r scheme is 
stable if (a) v = 0 and r > 0; or (b) v 2 = 1 and r > 1; or (c) 0 < v 2 < 1 and r = \v\ 
(Theorem 18). 

Excluding the four special cases addressed in Theorems 16 and 18, the set T defined 
in Eq. (4.66) is the set of all other (u, r) that satisfy the necessary stability conditions 
r > v 2 and v 1 < 1 (Theorem 19). To facilitate the development, T is divided into two 
disjoint subsets \k_ and which are defined in Eq. (4.66) and (4.67). 

It turns out that Eqs. (4.25) and (4.27) are satisfied by all (v. r) G T (Theorems 21 
and 22). Thus, according to Theorem 17, a given ( u , t ) G T is c-r stable if and only if 
it satisfies Eq. (4.26). As such, one arrives at the conclusion that a given (u, r) G df is 
c-r stable if and only if it satisfies Eq. (4.84) (Theorem 23). This necessary and sufficient 
stability condition obviously is even simpler than those given in Theorem 17. 

With the aid of Theorems 24-31, for the set T, we are able to obtain the explicit 
solution to the necessary and sufficient stability condition Eq. (4.84) in the form given in 
Theorem 32. The functions I+(x), J+(x), K + (x), L + (x), and L-(x), 0 < x < 1, that 
appear in Theorem 32 are defined in Eqs. (4.91), (4.92), (4.114), and (4.122). 

In principle, whether a given (za r) is c-r stable can be determined by using Theorems 
12, 16, 18, 19, and 32. However, by using the alternative definitions of and T+ given 
in Theorem 33, and the ordering properties Eqs. (4.160)-(4.168) given in Theorem 34, it 
was shown that Theorems 12, 16, 18, 19, and 32 can be combined and turned into the 
simple explicit form of necessary and sufficient stability conditions given in Theorem 35. 

Finally note that the proof of the ordering properties Eqs. (4.160)-(4.168) is hinged 
on the rather incredible facts that the 4-6th order polynomials in x or y/x that appear 
in Eqs. (4.247), (4.251), (4.255), (4.259), (4.263), (4.295), (4.303), (4.312), and (4.331) all 
can be factorized and studied analytically. 
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Appendix A. Numerical Validation of Theorem 34 
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implicit real*8 (a-h, o-z ) 
c 

c Program "ineqs.for". 

c 

c This program is used to verify numerically the inequalities 

c Equations ( 4 . 1 60 ) -- ( 4 . 1 68 ) (see Theorem 34). 

c 

c *** The functions I-plus, K-plus, L-minus, and L-plus are undefined 
c *** at x=0.d0. Thus, instead of being evaluated at x=0.d0, these 
c *** functions will be evaluated at x=ep where ep is a very small 
c *** positive number, 
c 

c *** At x=l.d0, 2 . dO* (2 . d0-x-x**2 ) =0 . dO . Because of round-off errors, 
c *** the value of this expression may become negative when x is very 
c *** close to 1 . dO . As such the square root of this expression and 
c *** therefore the functions L-minus and L-plus may be undefined 
c *** computationally when x is too close to 1 . dO . Thus, instead of 
c *** being evaluated at x=l.d0, the functions will be evaluated at 
c *** x=l.d0-eq where eq is a very small positive number, 
c 

c *** n l = number of uniform sub-intervals in ( 0 , cl) . 

c *** n2 = number of uniform sub-intervals in (cl,c2) . 

c *** n3 = number of uniform sub-intervals in (c2,c3) . 

c *** n 4 = number of uniform sub-intervals in (c3,c4) . 

c *** n5 = number of uniform sub-intervals in (c4,l) . 

c 

srt (x) =dsqrt (x) 

f ip (x) = (3.d0*x-2.d0+2. dO* dsqrt ( 3 . d0*x**2-3 . d0*x+l . dO ) ) /x 
f jp (x) = (3.d0*x-2 .dO+dsqrt (2 .dO* (x**3-x+2 . dO ) ) ) / (2 .dO-x) 
fkp (x) = (x— 1 .dO+dsqrt (5 . d0*x**2-2 . d0*x+l . dO ) ) / (2 .d0*x) 
flm(x) = (4 .dO-x-2 .d0*dsqrt (2 .dO* (2 . d0-x-x**2 ) ) ) /x 
flp (x) = (4 .dO-x+2 .d0*dsqrt (2 .dO* (2 . d0-x-x**2 ) ) ) /x 
c 

nl = 17 
n2 = 10 
n3=12 
n4 = 14 
n5=47 
ep=l .d-7 
eq=l . d-12 
one=l . dO-eq 
n5m=n5-l 
c 

cl=3 .dO-2 .d0*dsqrt (2 .dO) 
c2=3. dO/11. dO 

c3=(41.d0-7. dO+dsqrt ( 33 . dO ) ) /2 . dO 

c4= (dexp ( (1 .dO/3 .dO) ^dlog (dsqrt ( 1 664 . dO/27 . dO ) +181 . dO/27 . dO ) ) 

* -dexp ( (1 .dO/3 .dO) *dlog (dsqrt ( 1 664 . dO/27 . dO ) -181 . dO/27 . dO ) ) 

* -2 . dO/3 . dO ) ^^2 
c 

open (unit=8 , f ile= 1 ineqs . txt * ) 
write (8,1) 
write (8,2) 

write (8,3) nl , n2 , n3 , n4 , n5 
write (8,4) ep,eq 
write (8,2) 

write (8,10) cl,c2,c3,c4 
write (8,2) 
c 

dxl=cl/df loat (nl) 
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c 


100 


200 


300 


400 


500 


1 

2 

3 

4 

10 

20 

30 


dx2= ( c2-cl ) /dfloat (n2) 
dx3= (c3-c2) /dfloat (n3) 
dx4 = (c4-c3 ) /dfloat (n4) 
dx5= ( 1 . d0-c4 ) /dfloat (n5) 


write (8,20) ep 

write (8,30) f ip (ep) , ep, f kp (ep) , f lm (ep) 
write (8,40) f jp (ep) , srt (ep) , f lp (ep) 
write (8,2) 
x=0 . dO 

do 100 i=l,nl 
x=x+dxl 

write (8,20) x 

write (8,30) f ip (x) , x, f kp (x) , f lm (x) 

write (8,40) f jp (x) , srt (x) , f lp (x) 

continue 

write (8,2) 

x=cl 

do 200 i=l,n2 
x=x+dx2 

write (8,20) x 

write (8,50) x, f ip (x) , f kp (x) , f lm (x) 

write (8,40) f jp (x) , srt (x) , f lp (x) 

continue 

write (8,2) 

x=c2 

do 300 i=l,n3 
x=x+dx3 

write (8,20) x 

write (8,60) x, f kp (x) , f lm (x) , f ip (x) 

write (8,40) f jp (x) , srt (x) , f lp (x) 

continue 

write (8,2) 

x=c3 


do 400 i=l,n4 
x=x+dx4 


write (8,20) x 

write (8,70) x, f kp (x) , f lm (x) , srt (x) 

write (8,80) f jp (x) , f ip (x) , f lp (x) 

continue 

write (8,2) 

x=c4 

do 500 i=l,n5m 
x=x+dx5 


write (8,20) x 

write (8,90) x, f kp (x) , srt (x) , f lm (x) 
write (8,80) f jp (x) , f ip (x) , f lp (x) 
continue 
write (8,2) 

(8,20) one 

( 8 , 90 ) one, f kp (one) , srt (one) , f lm (one) 
(8,80) f jp (one) , f ip (one) , f lp (one) 
(unit=8 ) 


write 

write 

write 

close 

format 

format 

format 

format 

format 

format 

format 


***** ph e output for the code "ineqs.for". ***** 

nl = ’ , i3 , 1 n2 =',i3,' n3 =',i3,' n4 =',i3,' n5 = 
ep =' ,gl4.7, ' eq = ' ,gl4.7) 

cl =' ,gl4.7, ' c2 = 1 , gl4 . 7 , ' c3 = l ,gl4.7,’ c4 = 
x = ' ,gl4 .7) 

f ip = 1 , gl4 . 7 , ' x = 1 , gl4 . 7 , 1 fkp =',gl4.7,' flm 


****** 

, i3 ) 
gi4 . 7) 

' gl 4 . 7 ; 


) 

) 
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40 

format 

(' fj P = ' , gl4 . 7, 

' srt 

50 

format 

X 

II 

fip = 

60 

format 

(' x = ' , gl4 . 7, ' 

fkp = 

70 

format 

(' x = ' , gl4 . 7, ' 

fkp = 

80 

format 

(' fjp = ' , g!4 . 7, 

' fip 

90 

format 

stop 

end 

X 

II 

M 

4^ 

fkp = 


=',gl4.7. 

' fip 

=',gl4.7) 



,gl4.7, ' 

fkp = ' 

,gi4.7, 1 

flm 

= 1 gl4 , 

.7) 

, gl4 . 7, ' 

flm = ' 

, gl4 . 7, ' 

fip 

= 1 gl4 , 

■ 7) 

, gl4 . 7, ' 

flm = ' 

, gl4 . 7, ' 

srt 

= 1 gl4 , 

■ 7) 

=',gl4.7. 

' fip 

='gl4 .7) 




,gl4.7, ' 

srt = ' 

,gl4.7, ' 

flm 

= 1 gl4 , 

• 7) 
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x = 0.1009252E-01 

fip = 0.7685593E-02 x = 0 . 1009252E-01 fkp = 0. 1019436E-01 flm = 0 . 1138297E-01 

fjp = 0.1267669E-01 srt = 0.1004615 flp = 790.6547 

x = 0.2018504E-01 

fip = 0.1561019E-01 x = 0.2018504E-01 fkp = 0 . 2059214E-01 flm = 0 . 2282467E-01 

fjp = 0 . 2547771E-01 srt = 0.1420741 flp = 394.3102 

x = 0.3027757E-01 

fip = 0.2378408E-01 x = 0 . 3027757E-01 fkp = 0. 3119254E-01 flm = 0 . 3432655E-01 

fjp = 0.3840651E-01 srt = 0.1740045 flp = 262.1877 

x = 0.4037009E-01 

fip = 0.3221806E-01 x = 0.4037009E-01 fkp = 0 . 4199420E-01 flm = 0 . 4589012E-01 

fjp = 0.5146663E-01 srt = 0.2009231 flp = 196.1206 

x = 0.5046261E-01 

fip = 0.4092349E-01 x = 0 . 5046261E-01 fkp = 0. 5299516E-01 flm = 0 . 5751692E-01 

fjp = 0.6466170E-01 srt = 0.2246388 flp = 156.4757 

x = 0.6055513E-01 

fip = 0.4991230E-01 x = 0 . 6055513E-01 fkp = 0 . 6419281E-01 flm = 0 . 6920851E-01 

fjp = 0.7799544E-01 srt = 0.2460795 flp = 130.0418 

x = 0.7064765E-01 

fip = 0.5919705E-01 x = 0 . 7064765E-01 fkp = 0. 7558387E-01 flm = 0 . 8096653E-01 

fjp = 0.9147166E-01 srt = 0.2657963 flp = 111.1570 

x = 0.8074018E-01 

fip = 0.6879090E-01 x = 0 . 8074018E-01 fkp = 0 . 8716441E-01 flm = 0 . 9279262E-01 

fjp = 0.1050943 srt = 0.2841482 flp = 96.99047 

x = 0.9083270E-01 

fip = 0 . 7870770E-01 x = 0 . 9083270E-01 fkp = 0.9892977E-01 flm = 0.1046885 

fjp = 0.1188673 srt = 0.3013846 flp = 85.96932 

x = 0.1009252 

fip = 0.8896199E-01 x = 0.1009252 fkp = 0.1108746 flm = 0.1166559 

fjp = 0.1327948 srt = 0.3176873 flp = 77.14995 

x = 0.1110177 

fip = 0.9956905E-01 x = 0.1110177 fkp = 0.1229927 flm = 0.1286967 

fjp = 0.1468811 srt = 0.3331933 flp = 69.93186 

x = 0.1211103 

fip = 0.1105449 x = 0.1211103 fkp = 0.1352774 flm = 0.1408127 

fjp = 0.1611304 srt = 0.3480090 flp = 63.91469 

x = 0.1312028 

fip = 0.1219063 x = 0.1312028 fkp = 0.1477212 flm = 0.1530058 

fjp = 0.1755472 srt = 0.3622193 flp = 58.82131 

x = 0.1412953 

fip = 0.1336709 x = 0.1412953 fkp = 0.1603157 flm = 0.1652779 

fjp = 0.1901361 srt = 0.3758927 flp = 54.45373 

x = 0.1513878 

fip = 0.1458572 x = 0.1513878 fkp = 0.1730522 flm = 0.1776312 

fjp = 0.2049017 srt = 0.3890859 flp = 50.66677 

x = 0.1614804 

fip = 0.1584845 x = 0.1614804 fkp = 0.1859211 flm = 0.1900676 

fjp = 0.2198489 srt = 0.4018462 flp = 47.35156 

x = 0.1715729 


fip = 0.1715729 
fjp = 0.2349824 

‘Jr *Jr »Jr *Jr 


x = 0.1715729 
srt = 0.4142136 

r'r *■'** *•'** r'r 


fkp = 0.1989124 
flp = 44.42483 

*•*.- **'- **.V **'- **.V ‘Jr *Jr *Jr *Jr *Jr *Jr 


flm = 0.2025894 
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i neqs . txt 


x = 0. 1816883 
x = 0.1816883 
fjp = 0.2503425 
x = 0.1918038 
x = 0.1918038 
fjp = 0.2659001 
x = 0.2019192 
x = 0.2019192 
fjp = 0.2816607 
x = 0.2120346 
x = 0.2120346 
fjp = 0.2976296 
x = 0.2221501 
x = 0.2221501 
fjp = 0.3138126 
x = 0.2322655 
x = 0.2322655 
fjp = 0.3302153 
x = 0.2423810 
x = 0.2423810 
fjp = 0.3468438 
x = 0.2524964 
x = 0.2524964 
fjp = 0.3637040 
x = 0.2626118 
x = 0.2626118 
fjp = 0.3808023 
x = 0.2727273 
x = 0.2727273 
fjp = 0.3981449 


x = 0.2828359 
x = 0.2828359 
fjp = 0.4157266 
x = 0.2929445 
x = 0.2929445 
fjp = 0.4335657 
x = 0.3030531 
x = 0.3030531 
fjp = 0.4516690 
x = 0.3131618 
x = 0.3131618 
fjp = 0.4700437 
x = 0.3232704 
x = 0.3232704 
fjp = 0.4886969 
x = 0.3333790 
x = 0.3333790 
fjp = 0.5076361 
x = 0.3434876 
x = 0.3434876 
fjp = 0.5268689 
x = 0.3535962 
x = 0.3535962 
fjp = 0.5464032 
x = 0.3637049 
x = 0.3637049 
fjp = 0.5662470 
x = 0.3738135 
x = 0.3738135 
fjp = 0.5864085 
x = 0.3839221 
x = 0.3839221 


fip = 0.1851749 
srt = 0.4262491 

fip = 0.1992835 
srt = 0.4379541 

fip = 0.2139216 
srt = 0.4493542 

fip = 0.2291133 
srt = 0.4604722 

fip = 0.2448834 
srt = 0.4713280 

fip = 0.2612574 
srt = 0.4819393 

fip = 0.2782619 
srt = 0.4923220 

fip = 0.2959241 
srt = 0.5024902 

fip = 0.3142718 
srt = 0.5124567 

fip = 0.3333333 
srt = 0.5222330 


fkp = 0.3469168 
srt = 0.5318232 

fkp = 0.3604781 
srt = 0.5412435 

fkp = 0.3740056 
srt = 0.5505026 

fkp = 0.3874879 
srt = 0.5596086 

fkp = 0.4009141 
srt = 0.5685687 

fkp = 0.4142738 
srt = 0.5773898 

fkp = 0.4275569 
srt = 0.5860782 

fkp = 0.4407542 
srt = 0.5946396 

fkp = 0.4538567 
srt = 0.6030795 

fkp = 0.4668562 
srt = 0.6114029 

fkp = 0.4797451 


fkp = 0.2120452 
fip = 41.81622 

fkp = 0.2252789 
fip = 39.48134 

fkp = 0.2386021 
fip = 37.37903 

fkp = 0.2520026 
fip = 35.47599 

fkp = 0.2654682 
fip = 33.74499 

fkp = 0.2789864 
fip = 32.16352 

fkp = 0.2925447 
fip = 30.71286 

fkp = 0.3061303 
fip = 29.37726 

fkp = 0.3197307 
fip = 28.14342 

fkp = 0.3333333 
fip = 27.00000 


flm = 0.3469816 
fip = 25.93797 

flm = 0.3607478 
fip = 24.94818 

flm = 0.3746351 
fip = 24.02338 

flm = 0.3886471 
fip = 23.15726 

flm = 0.4027873 
fip = 22.34430 

flm = 0.4170595 
fip = 21.57965 

flm = 0.4314677 
fip = 20.85904 

flm = 0.4460158 
fip = 20.17866 

flm = 0.4607081 
fip = 19.53515 

flm = 0.4755489 
fip = 18.92550 

flm = 0.4905430 


flm = 0.2152275 
flm = 0.2279558 
flm = 0.2407767 
flm = 0.2536927 
flm = 0.2667063 
flm = 0.2798201 
flm = 0.2930369 
flm = 0.3063594 
flm = 0.3197905 
flm = 0.3333333 

fip = 0.3531240 
fip = 0.3736852 
fip = 0.3950460 
fip = 0.4172353 
fip = 0.4402814 
fip = 0.4642118 
fip = 0.4890530 
fip = 0.5148299 
fip = 0.5415657 
fip = 0.5692811 
fip = 0.5979939 


* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
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i neqs.txt 


fjp 

= 0.6068962 

srt 

= 0.6196145 

fip 

= 18.34702 




X = 

0.3940307 








X = 

0.3940307 

fkp = 

0.4925164 

flm = 

0.5056950 

fip 

= 

0.6277187 

fjp, 

= 0.6277187 

srt 

=0.6277187 

..fip 

= 17.79729 




X = 

0.4037587 








X = 

0.4037587 

fkp = 

0.5046894 

flm = 

0.5204302 

srt 

= 

0.6354201 

fjp 

= 0.6480814 

fip 

= 0.6572896 

fip 

= 17.29339 




X = 

0.4134866 








X = 

0.4134866 

fkp = 

0.5167423 

flm = 

0.5353209 

srt 

= 

0.6430292 

fjp 

= 0.6687705 

fip 

= 0.6878137 

fip 

= 16.81234 




X = 

0.4232145 








X = 

0.4232145 

fkp = 

0.5286703 

flm = 

0.5503721 

srt 

= 

0.6505494 

fjp 

= 0.6897943 

fip 

= 0.7192926 

fip 

= 16.35257 




X = 

0.4329424 








X = 

0.4329424 

fkp = 

0.5404690 

flm = 

0.5655889 

srt 

= 

0.6579836 

fjp 

= 0.7111613 

fip 

= 0.7517238 

fip 

= 15.91262 




X = 

0.4426703 








X = 

0.4426703 

fkp = 

0.5521347 

flm = 

0.5809764 

srt 

= 

0.6653348 

fjp 

= 0.7328803 

fip 

= 0.7851000 

fip 

= 15.49116 




X = 

0.4523983 








X = 

0.4523983 

fkp = 

0.5636642 

flm = 

0.5965403 

srt 

= 

0.6726056 

fjp 

= 0.7549601 

fip 

= 0.8194089 

fip 

= 15.08699 




X = 

0.4621262 








X = 

0.4621262 

fkp = 

0.5750545 

flm = 

0.6122865 

srt 

= 

0.6797986 

fjp 

= 0.7774099 

fip 

= 0.8546332 

fip 

= 14.69900 




X = 

0.4718541 








X = 

0.4718541 

fkp = 

0.5863032 

flm = 

0.6282209 

srt 

= 

0.6869164 

fjp 

= 0.8002391 

fip 

= 0.8907504 

fip 

= 14.32617 




X = 

0.4815820 








X = 

0.4815820 

fkp = 

0.5974081 

flm = 

0.6443499 

srt 

= 

0.6939611 

fjp 

= 0.8234573 

fip 

= 0.9277326 

fip 

= 13.96757 




X = 

0.4913100 








X = 

0.4913100 

fkp = 

0.6083677 

flm = 

0.6606804 

srt 

= 

0.7009351 

fjp 

= 0.8470743 

fip 

= 0.9655470 

fip 

= 13.62232 




X = 

0.5010379 








X = 

0.5010379 

fkp = 

0.6191806 

flm = 

0.6772194 

srt 

= 

0.7078403 

fjp 

= 0.8711002 

fip 

= 1.004156 

fip 

= 13.28964 




X = 

0.5107658 








X = 

0.5107658 

fkp = 

0.6298458 

flm = 

0.6939742 

srt 

= 

0.7146788 

fjp 

= 0.8955453 

fip 

= 1.043517 

fip 

= 12.96878 




X = 

0.5204937 








X = 

0.5204937 

fkp = 

0.6403626 

flm = 

0.7109527 

srt 

= 

0.7214525 

fjp 

= 0.9204201 

fip 

= 1.083583 

fip 

= 12.65907 




X = 

0.5302216 








X = 

0.5302216 

fkp = 

0.6507306 

flm = 

0.7281632 

srt 

= 

0.7281632 

fjp 

= 0.9457355 

fip 

= 1.124304 

fip 

= 12.35987 




X = 

0.5402169 








X = 

0.5402169 

fkp = 

0.6612284 

srt = 

0.7349945 

flm 

= 

0.7460974 

fjp 

= 0.9722173 

fip 

= 1.166769 

fip 

= 12.06277 




X = 

0.5502122 








X = 

0.5502122 

fkp = 

0.6715692 

srt = 

0.7417629 

flm 

= 

0.7642956 

fjp 

= 0.9991881 

fip 

= 1.209809 

fip 

= 11.77555 




X = 

0.5602075 








X = 

0.5602075 

fkp = 

0.6817533 

srt = 

0.7484701 

flm 

= 

0.7827683 

fjp 

= 1.026661 

fip 

= 1.253358 

fip 

= 11.49765 




X = 

0.5702028 








X = 

0.5702028 

fkp = 

0.6917814 

srt = 

0.7551177 

flm 

= 

0.8015269 

fjp 

= 1.054648 

fip 

= 1.297353 

fip 

= 11.22857 




X = 

0.5801981 








X = 

0.5801981 

fkp = 

0.7016542 

srt = 

0.7617073 

flm 

= 

0.8205830 

fjp 

= 1.083164 

fip 

= 1.341727 

fip 

= 10.96781 
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i neqs . txt 


x = 0.5901934 
x = 0.5901934 
fjp = 1.112220 

x = 0.6001886 
x = 0.6001886 
fjp = 1.141833 

x = 0.6101839 
x = 0.6101839 
fjp = 1.172015 

x = 0.6201792 
x = 0.6201792 
fjp = 1.202782 

x = 0.6301745 
x = 0.6301745 
fjp = 1.234148 

x = 0.6401698 
x = 0.6401698 
fjp = 1.266130 

x = 0.6501651 
x = 0.6501651 
fjp = 1.298742 

x = 0.6601603 
x = 0.6601603 
fjp = 1.332001 

x = 0.6701556 
x = 0.6701556 
fjp = 1.365924 

x = 0.6801509 
x = 0.6801509 
fjp = 1.400528 

x = 0.6901462 
x = 0.6901462 
fjp = 1.435830 

x = 0.7001415 
x = 0.7001415 
fjp = 1.471848 

x = 0.7101368 
x = 0.7101368 
fjp = 1.508601 

x = 0.7201320 
x = 0.7201320 
fjp = 1.546108 

x = 0.7301273 
x = 0.7301273 
fjp = 1.584388 

x = 0.7401226 
x = 0.7401226 
fjp = 1.623462 

x = 0.7501179 
x = 0.7501179 
fjp = 1.663349 

x = 0.7601132 
x = 0.7601132 
fjp = 1.704072 

x = 0.7701085 
x = 0.7701085 
fjp = 1.745652 

x = 0.7801038 
x = 0.7801038 
fjp = 1.788111 

x = 0.7900990 
x = 0.7900990 
fjp = 1.831472 


fkp = 0.7113728 
fip = 1.386415 

fkp = 0.7209383 
fip = 1.431351 

fkp = 0.7303522 
fip = 1.476471 

fkp = 0.7396159 
fip = 1.521713 

fkp = 0.7487310 
fip = 1.567017 

fkp = 0.7576992 
fip = 1.612325 

fkp = 0.7665224 
fip = 1.657582 

fkp = 0.7752023 
fip = 1.702736 

fkp = 0.7837410 
fip = 1.747740 

fkp = 0.7921404 
fip = 1.792548 

fkp = 0.8004027 
fip = 1.837119 

fkp = 0.8085299 
fip = 1.881414 

fkp = 0.8165241 
fip = 1.925400 

fkp = 0.8243875 
fip = 1.969045 

fkp = 0.8321223 
fip = 2.012320 

fkp = 0.8397306 
fip = 2.055201 

fkp = 0.8472146 
fip = 2.097666 

fkp = 0.8545766 
fip = 2.139695 

fkp = 0.8618186 
fip = 2.181271 

fkp = 0.8689429 
fip = 2.222381 

fkp = 0.8759515 
fip = 2.263011 


srt = 0.7682404 
fip = 10.71493 

srt = 0.7747184 
fip = 10.46950 

srt = 0.7811427 
fip = 10.23113 

srt = 0.7875146 
fip = 9.999448 

srt = 0.7938353 
fip = 9.774095 

srt = 0.8001061 
fip = 9.554745 

srt = 0.8063281 
fip = 9.341082 

srt = 0.8125025 
fip = 9.132810 

srt = 0.8186303 
fip = 8.929647 

srt = 0.8247126 
fip = 8.731324 

srt = 0.8307504 
fip = 8.537585 

srt = 0.8367446 
fip = 8.348183 

srt = 0.8426961 
fip = 8.162884 

srt = 0.8486059 
fip = 7.981461 

srt = 0.8544749 
fip = 7.803693 

srt = 0.8603038 
fip = 7.629368 

srt = 0.8660935 
fip = 7.458277 

srt = 0.8718447 
fip = 7.290217 

srt = 0.8775582 
fip = 7.124987 

srt = 0.8832348 
fip = 6.962386 

srt = 0.8888752 
fip = 6.802214 


flm = 0.8399495 
flm = 0.8596397 
flm = 0.8796680 
flm = 0.9000497 
flm = 0.9208013 
flm = 0.9419404 
flm = 0.9634858 
flm = 0.9854579 
flm = 1.007879 

flm = 1.030771 

flm = 1.054162 

flm = 1.078079 

flm = 1.102551 

flm = 1.127613 

flm = 1.153300 

flm = 1.179652 

flm = 1.206713 

flm = 1.234531 

flm = 1.263160 

flm = 1.292660 

flm = 1.323099 
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x = 0.8000943 


X = 

0.8000943 

fkp = 

0.8828468 

srt = 

0.8944799 

flm = 

1.354550 

fjp 

= 1.875760 

fip 

= 2.303152 

fip 

= 6.644271 



X = 

0.8100896 





X = 

0.8100896 

fkp = 

0.8896306 

srt = 

0.9000498 

flm = 

1.387102 

fjp 

= 1.920998 

fip 

= 2.342796 

fip 

= 6.488349 



X = 

0.8200849 





X = 

0.8200849 

fkp = 

0.8963053 

srt = 

0.9055854 

flm = 

1.420850 

fjp 

= 1.967213 

fip 

= 2.381936 

fip 

= 6.334238 



X = 

0.8300802 





X = 

0.8300802 

fkp = 

0.9028728 

srt = 

0.9110874 

flm = 

1.455906 

fjp 

= 2.014428 

fip 

= 2.420567 

fip 

= 6.181717 



X = 

0.8400755 





X = 

0.8400755 

fkp = 

0.9093352 

srt = 

0.9165563 

flm = 

1.492400 

fjp 

= 2.062673 

fip 

= 2.458687 

fip 

= 6.030554 



X = 

0.8500707 





X = 

0.8500707 

fkp = 

0.9156945 

srt = 

0.9219928 

flm = 

1.530482 

fjp 

= 2.111973 

fip 

= 2.496292 

fip 

= 5.880499 



X = 

0.8600660 





X = 

0.8600660 

fkp = 

0.9219528 

srt = 

0.9273974 

flm = 

1.570329 

fjp 

= 2.162358 

fip 

= 2.533382 

fip 

= 5.731282 



X = 

0.8700613 





X = 

0.8700613 

fkp = 

0.9281119 

srt = 

0.9327708 

flm = 

1.612151 

fjp 

= 2.213857 

fip 

= 2.569957 

fip 

= 5.582603 



X = 

0.8800566 





X = 

0.8800566 

fkp = 

0.9341739 

srt = 

0.9381133 

flm = 

1.656201 

fjp 

= 2.266500 

fip 

= 2.606019 

fip 

= 5.434124 



X = 

0.8900519 





X = 

0.8900519 

fkp = 

0.9401406 

srt = 

0.9434256 

flm = 

1.702787 

fjp 

= 2.320318 

fip 

= 2.641569 

fip 

= 5.285453 



X = 

0.9000472 





X = 

0.9000472 

fkp = 

0.9460140 

srt = 

0.9487082 

flm = 

1.752292 

fjp 

= 2.375344 

fip 

= 2.676609 

fip 

= 5.136132 



X = 

0.9100424 





X = 

0.9100424 

fkp = 

0.9517957 

srt = 

0.9539614 

flm = 

1.805199 

fjp 

= 2.431612 

fip 

= 2.711144 

fip 

= 4.985600 



X = 

0.9200377 





X = 

0.9200377 

fkp = 

0.9574878 

srt = 

0.9591860 

flm = 

1.862136 

fjp 

= 2.489155 

fip 

= 2.745178 

fip 

= 4.833160 



X = 

0.9300330 





X = 

0.9300330 

fkp = 

0.9630919 

srt = 

0.9643822 

flm = 

1.923937 

fjp 

= 2.548009 

fip 

= 2.778714 

fip 

= 4.677908 



X = 

0.9400283 





X = 

0.9400283 

fkp = 

0.9686098 

srt = 

0.9695506 

flm = 

1.991756 

fjp 

= 2.608213 

fip 

= 2.811757 

fip 

= 4.518626 



X = 

0.9500236 





X = 

0.9500236 

fkp = 

0.9740431 

srt = 

0.9746915 

flm = 

2.067265 

fjp 

= 2.669803 

fip 

= 2.844314 

fip 

= 4.353579 



X = 

0.9600189 





X = 

0.9600189 

fkp = 

0.9793937 

srt = 

0.9798055 

flm = 

2.153047 

fjp 

= 2.732819 

fip 

= 2.876389 

fip 

= 4.180123 



X = 

0.9700141 





X = 

0.9700141 

fkp = 

0.9846630 

srt = 

0.9848930 

flm = 

2.253481 

fjp 

= 2.797304 

fip 

= 2.907989 

fip 

= 3.993822 



X = 

0.9800094 





X = 

0.9800094 

fkp = 

0.9898528 

srt = 

0.9899543 

flm = 

2.377166 

fjp 

= 2.863299 

fip 

= 2.939120 

fip 

= 3.786021 



X = 

0.9900047 





X = 

0.9900047 

fkp = 

0.9949646 

srt = 

0.9949898 

flm = 

2.546482 

fjp, 

= 2.930849 

fip. 

= 2.969788 

, ...fip, 

= 3.534287 



n « n i 

X = 
X = 

» « « « « « « « « n <V « « < 

1.000000 

1.000000 

fkp = 

* « « « « « « « « « « « i 

1.000000 

* « « « « n i 

srt = 

» « n « « « n « « « « 

1.000000 

flm = 

2.999995 
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fjp = 3.000000 fip = 3.000000 flp = 


3.000005 
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Appendix B. Numerical Validation of Theorem 35 
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implicit real*8 (a-h, o-z ) 

complex* 1 6 a, b, c, cdsqrt , xl , x2 , x3 , dcmplx 
c 

c Program "ctausc . f or " . 

c 

c This program is used to verify numerically the assertion made 

c in part G of Theorem 35. 

c 

c *** The critical value of tau, i . e . , tauo(nu**2), is evaluated for 
c *** each given value of nu (the Courant number) . 
c 

c *** Given any (nu^ tau) , the spectral radius of the amplification 
c *** matrix is a function of the phase angle theta. The least upper 
c *** bound (denoted by "am") of the spectral radii over the range 
c *** -pi < theta . le . pi is evaluated for each given (nu, tau) . 
c 

c *** When nu is replaced by -nu, each of the two resulting 
amplification 

c *** factors (defined in Eq. (4.7)) becomes the complex conjugate of 
that 

c *** before sign-change. Thus the spectral radius does not change as nu 
c *** is replaced by -nu . For this reason, the range of nu can be 
limited 

c *** to nu . ge . 0 . 
c 

c *** When theta is replaced by -theta, each of the two resulting 
c *** amplification factors also becomes the complex conjugate of 
c *** that before sign-change. Thus the range of theta can be limited 
c *** to 0 .le. theta . le . pi. 
c 

c *** Theorems 16 and 18 imply that the least upper bound am = 1 if 

c *** nu = 1 and tau . ge . 1 (Note: According to Eq. (4.7), the value 

c *** of the principal amplification factor = 1 when theta = 0. Thus 

c *** am . ge . 1 for any (nu, tau) . In turn, this implies that am = 1 

c *** for any (nu, tau) which meets the condition Eq. (4.2)) . Moreover, 
c *** Theorems 6 and 12 imply that am > 1 if nu > 1 regardless the 
c *** value assumed by tau. Thus numerical results may not be consistent 
c *** with theoretical predictions at the singular case nu = 1 if 
c *** round-off errors are not controlled carefully. For this reason, 
c *** a statement "if (dabs (x-1 . dO ) . lt.ep) x=l.d0" is added in the code 
c *** to insure that the value of x is really "1" as intended. Here 
c *** ep (>0) is an input parameter and assumes to be very small, 
c 

c x = nu . 

c z = The phase angle theta of a Fourier component, 

c nx = number of the values of nu . 

c 

c nt = number of the values of tau with tau>tauo (tau<tauo) for 

c each value of nu . Here tauo is the critical value of tau 

c associated with a given value of nu . Because the case with 

c tau=tauo is always considered, there are (2*nt+l) values 

c of tau associated with each value of nu, i.e., 

c tauo* (l-dt*nt) , tauo* ( 1-dt* (nt-1 )),... , tauo* (1-dt) , tauo, 

c tauo* (1+dt) , ... , tauo* (1+dt* (nt-1) ) , tauo* (l+dt*n) . 

c 

c nz = number of the intervals over the domain 

c 0 .le. theta . le . pi. 

c xs = The initial value of nu . 

c 

fkp (s) = ( s — 1 . dO+dsqrt (5 . d0*s**2-2 . d0*s + l . dO ) ) / (2 .d0*s) 
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c 


100 

200 

10 

15 


flm(s) = (4 . dO-s-2 . d0*dsqrt (2 . dO* (2.d0-s-s**2) ) ) /s 

pi = 3. 1415926535897932 dO 

nx = 2 5 

nt = 5 

nz = 1000 

xs = 0 . dO 

dx = 5 . d-2 

dt = i.d-4 

ep = l.d-7 

c2 = 3. dO/11. dO 

dz = pi/df loat (nz ) 

x = xs-dx 

nzp = nz+1 

ts = 1 . d0-dt*df loat (nt+1 ) 
nt2p = nt*2+l 

open (unit=8 , f ile- 1 ctausc . txt 1 ) 

write ( 8 , 10) 

write ( 8 , 15 ) 

write (8,20) nx,nt,nz 

write (8,30) xs,dx,dt,ep 

write ( 8 , 15 ) 

do 200 i = l,nx 

x = x+dx 

if (dabs (x-1 . dO ) . It . ep) x=l.d0 
xx = x**2 

if (xx.eq.O.dO) tauo = 0 . dO 

if (xx . gt . 0 . dO . and . xx . le . c2 ) tauo = flm(xx) 

if (xx.gt.c2) tauo = f kp (xx) 

tau = tauo*ts 

dtau = tauo*dt 

do 200 j = l,nt2p 

tau = tau+dtau 

am = 0 . dO 

z = -dz 

do 100 k = l,nzp 
z = z+dz 

zl = dcos ( z/2 . dO ) 
z2 = dsin ( z/2 . dO ) 
ar = l.dO+tau 
ai = 0 . dO 
br = -2.d0*tau*zl 
bi = x* (3.d0+tau) *z2 

cr = - ( ( 1 . dO - tau)*zl**2 + (l.dO + x**2)*z2**2) 
ci = -x*(l.d0 + tau) *zl*z2 
a = dcmplx (ar , ai) 
b = dcmplx (br, bi ) 
c = dcmplx (cr, ci) 

xl = (-b + cdsqrt (b**2 - 4 . d0*a*c) ) / (2 . d0*a) 

x2 = (-b - cdsqrt (b**2 - 4 . d0*a*c) ) / (2 . d0*a) 

al = cdabs (xl ) 

a2 = cdabs (x2 ) 

am = dmaxl (al , a2 , am) 

continue 

write (8,40) x, tauo, tau, am 

continue 

close (unit=8) 

format (' ***** The output for the code "ctausc . for" . ******') 
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20 

format 

c 

nx = ' , i4 , ' nt 

=',14,' nz = ' , i4 ) 


30 

format 

c 

xs ,gl4 .7, ' 

dx = ' ,gl4.7, ' dt = ' , 

gl4 . 7, ' ep 

40 

format 

c 

nu = ' , gll . 4 , ' 

tauo = f ,gl4.7,' tau 

=' ,gl4.7. 


* 

! 

am =' , g21 . 14 ) 




stop 

end 
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* * * * * 


ctauscBText . txt 

The output for the code "ctausc.for" . ****** 
nx = 25 nt = 5 nz =1000 

xs = 0.000000 dx = 0.5000000E-01 dt = 0.1000000E-03 ep 

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 


nu 

= 

0.000 

tauo 

nu 

= 

0.000 

tauo 

nu 

= 

0.000 

tauo 

nu 

= 

0.000 

tauo 

nu 

= 

0.000 

tauo 

nu 

= 

0.000 

tauo 

nu 

= 

0.000 

tauo 

nu 

- 

0.000 

tauo 

nu 

= 

0.000 

tauo 

nu 

= 

0.000 

tauo 

nu 

= 

0.000 

tauo 

nu 

= 

0.5000E-01 

tauo 

nu 

= 

0.5000E-01 

tauo 

nu 

= 

0.5000E-01 

tauo 

nu 

= 

0.5000E-01 

tauo 

nu 

= 

0.5000E-01 

tauo 

nu 

= 

0.5000E-01 

tauo 

nu 

= 

0.5000E-01 

tauo 

nu 

= 

0.5000E-01 

tauo 

nu 

= 

0.5000E-01 

tauo 

nu 

= 

0.5000E-01 

tauo 

nu 

= 

0.5000E-01 

tauo 

nu 

= 

0 . 1000 

tauo 

nu 

= 

0.1000 

tauo 

nu 

= 

0 . 1000 

tauo 

nu 

= 

0.1000 

tauo 

nu 

= 

0 . 1000 

tauo 

nu 

= 

0.1000 

tauo 

nu 

= 

0 . 1000 

tauo 

nu 

= 

0.1000 

tauo 

nu 

= 

0 . 1000 

tauo 

nu 

= 

0.1000 

tauo 

nu 

= 

0 . 1000 

tauo 

nu 

= 

0.1500 

tauo 

nu 

= 

0.1500 

tauo 

nu 

= 

0.1500 

tauo 

nu 

= 

0.1500 

tauo 

nu 

= 

0.1500 

tauo 

nu 

= 

0.1500 

tauo 

nu 

= 

0.1500 

tauo 

nu 

= 

0.1500 

tauo 

nu 

= 

0.1500 

tauo 

nu 

= 

0.1500 

tauo 

nu 

= 

0.1500 

tauo 

nu 

= 

0.2000 

tauo 

nu 

= 

0.2000 

tauo 

nu 

= 

0.2000 

tauo 

nu 

= 

0.2000 

tauo 

nu 

= 

0.2000 

tauo 

nu 

= 

0.2000 

tauo 

nu 

= 

0.2000 

tauo 

nu 

= 

0.2000 

tauo 

nu 

= 

0.2000 

tauo 

nu 

= 

0.2000 

tauo 

nu 

= 

0.2000 

tauo 

nu 

= 

0.2500 

tauo 

nu 

= 

0.2500 

tauo 

nu 

= 

0.2500 

tauo 


= 0 . 000000 tau 
= 0 . 000000 tau 
= 0 . 000000 tau 
= 0 . 000000 tau 
= 0 . 000000 tau 
= 0 . 000000 tau 
= 0 . 000000 tau 
= 0 . 000000 tau 
= 0 . 000000 tau 
= 0 . 000000 tau 
= 0 . 000000 tau 
= 0.2814261E-02 tau 
= 0.2814261E-02 tau 
= 0.2814261E-02 tau 
= 0.2814261E-02 tau 
= 0.2814261E-02 tau 
= 0.2814261E-02 tau 
= 0.2814261E-02 tau 
= 0.2814261E-02 tau 
= 0.2814261E-02 tau 
= 0.2814261E-02 tau 
= 0.2814261E-02 tau 
= 0.1127835E-01 tau 
= 0.1127835E-01 tau 
= 0.1127835E-01 tau 
= 0.1127835E-01 tau 
= 0.1127835E-01 tau 
= 0.1127835E-01 tau 
= 0.1127835E-01 tau 
= 0.1127835E-01 tau 
= 0.1127835E-01 tau 
= 0.1127835E-01 tau 
= 0.1127835E-01 tau 
= 0.2545752E-01 tau 
= 0.2545752E-01 tau 
= 0.2545752E-01 tau 
= 0.2545752E-01 tau 
= 0.2545752E-01 tau 
= 0.2545752E-01 tau 
= 0.2545752E-01 tau 
= 0.2545752E-01 tau 
= 0.2545752E-01 tau 
= 0.2545752E-01 tau 
= 0.2545752E-01 tau 
= 0.4546499E-01 tau 
= 0.4546499E-01 tau 
= 0.4546499E-01 tau 
= 0.4546499E-01 tau 
= 0.4546499E-01 tau 
= 0.4546499E-01 tau 
= 0.4546499E-01 tau 
= 0.4546499E-01 tau 
= 0.4546499E-01 tau 
= 0.4546499E-01 tau 
= 0.4546499E-01 tau 
= 0.7146911E-01 tau 
= 0.7146911E-01 tau 
= 0.7146911E-01 tau 


= 0.000000 am 
= 0.000000 am 
= 0.000000 am 
= 0.000000 am 
= 0.000000 am 
= 0.000000 am 
= 0.000000 am 
= 0.000000 am 
= 0.000000 am 
= 0.000000 am 
= 0.000000 am 
= 0.2812854E-02 am 
= 0.2813136E-02 am 
= 0.2813417E-02 am 
= 0.2813699E-02 am 
= 0.2813980E-02 am 
= 0.2814261E-02 am 
= 0.2814543E-02 am 
= 0.2814824E-02 am 
= 0.2815106E-02 am 
= 0.2815387E-02 am 
= 0.2815669E-02 am 
= 0.1127272E-01 am 
= 0.1127384E-01 am 
= 0.1127497E-01 am 
= 0.1127610E-01 am 
= 0.1127723E-01 am 
= 0.1127835E-01 am 
= 0.1127948E-01 am 
= 0.1128061E-01 am 
= 0.1128174E-01 am 
= 0. 1128287E-01 am 
= 0.1128399E-01 am 
= 0.2544479E-01 am 
= 0.2544734E-01 am 
= 0.2544988E-01 am 
= 0.2545243E-01 am 
= 0.2545498E-01 am 
= 0.2545752E-01 am 
= 0 . 2546007E-01 am 
= 0.2546261E-01 am 
= 0.2546516E-01 am 
= 0.2546770E-01 am 
= 0.2547025E-01 am 
= 0.4544225E-01 am 
= 0.4544680E-01 am 
= 0.4545135E-01 am 
= 0.4545589E-01 am 
= 0.4546044E-01 am 
= 0.4546499E-01 am 
= 0.4546953E-01 am 
= 0.4547408E-01 am 
= 0.4547862E-01 am 
= 0.4548317E-01 am 
= 0.4548772E-01 am 
= 0.7143338E-01 am 
= 0.7144052E-01 am 
= 0.7144767E-01 am 
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nu 

= 

0.2500 

tauo 

= 

0.7146911E-01 

tau 

- 

0.7145482E-01 

am 

= 

1.0000011789100 

nu 

= 

0.2500 

tauo 

= 

0.7146911E-01 

tau 

= 

0.7146196E-01 

am 

= 

1.0000005880919 

nu 

= 

0.2500 

tauo 

= 

0.7146911E-01 

tau 

= 

0.7146911E-01 
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= 

1.0000000000000 

nu 
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0.2500 

tauo 

= 

0.7146911E-01 

tau 

= 

0.7147626E-01 
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1.0000000000000 

nu 

= 

0.2500 

tauo 

= 

0.7146911E-01 

tau 

= 

0.7148340E-01 
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1.0000000000000 

nu 

= 

0.2500 

tauo 

= 

0.7146911E-01 

tau 

= 

0.7149055E-01 
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= 

1.0000000000000 

nu 

= 

0.2500 

tauo 

= 

0.7146911E-01 

tau 

= 

0.7149770E-01 
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= 

1.0000000000000 

nu 

= 

0.2500 

tauo 

= 

0.7146911E-01 

tau 

= 

0.7150484E-01 
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= 

1.0000000000000 

nu 

= 

0.3000 

tauo 

= 

0.1037043 

tau 

= 

0.1036525 
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= 

1.0000034247159 

nu 

= 

0.3000 

tauo 

= 

0.1037043 

tau 

= 

0.1036629 
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= 

1.0000027345694 

nu 

= 

0.3000 

tauo 

= 

0.1037043 

tau 

= 

0.1036732 
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= 

1.0000020464043 

nu 

= 

0.3000 

tauo 

= 

0.1037043 

tau 

= 

0.1036836 

am 

= 

1.0000013620442 

nu 

= 

0.3000 

tauo 

= 

0.1037043 

tau 

= 

0.1036940 
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= 

1.0000006790801 

nu 

= 

0.3000 

tauo 

= 

0.1037043 
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= 

0.1037043 
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= 

1.0000000000000 

nu 

= 

0.3000 

tauo 

= 

0.1037043 

tau 

= 

0.1037147 
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= 

1.0000000000000 

nu 

= 

0.3000 

tauo 

= 

0.1037043 

tau 

= 

0.1037251 

am 

= 

1.0000000000000 

nu 

= 

0.3000 

tauo 

= 

0.1037043 

tau 

= 

0.1037354 
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1.0000000000000 

nu 

= 

0.3000 

tauo 

= 

0.1037043 

tau 

= 

0.1037458 

am 

= 

1.0000000000000 

nu 

= 

0.3000 

tauo 

= 

0.1037043 

tau 

= 

0.1037562 

am 

= 

1.0000000000000 

nu 

= 

0.3500 

tauo 

= 

0.1424870 
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= 

0.1424158 

am 

= 

1.0000034155413 

nu 
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0.3500 

tauo 

= 

0.1424870 
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0.1424301 
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1.0000027246816 
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0.3500 
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0.1424870 
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1.0000020381007 
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= 

0.3500 
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= 

0.1424870 
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= 

0.1424586 
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= 

1.0000013549685 

nu 

= 

0.3500 

tauo 

= 

0.1424870 

tau 

= 

0.1424728 

am 

= 

1.0000006753657 

nu 

= 

0.3500 

tauo 

= 

0.1424870 

tau 

= 

0.1424870 

am 

= 

1.0000000000000 

nu 

= 

0.3500 

tauo 

= 

0.1424870 

tau 

= 

0.1425013 

am 

= 

1.0000000000000 

nu 

= 

0.3500 

tauo 

= 

0.1424870 

tau 

= 

0.1425155 

am 

= 

1.0000000000000 

nu 

= 

0.3500 

tauo 

= 

0.1424870 

tau 

= 

0.1425298 

am 

= 

1.0000000000000 

nu 

= 

0.3500 

tauo 

= 

0.1424870 

tau 

= 

0.1425440 

am 

= 

1.0000000000000 

nu 

= 

0.3500 

tauo 

= 

0.1424870 

tau 

= 

0.1425583 

am 

= 

1.0000000000000 
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= 

0.4000 
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= 

0.1882382 

tau 

= 

0.1881441 

am 

= 

1.0000027984482 

nu 

= 

0.4000 
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= 

0.1882382 
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= 
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1.0000022284014 
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0.4000 
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= 

0.1882382 
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= 

0.4000 
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= 

0.1882382 
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= 
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tauo 
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0.1882382 
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0.4500 
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0.2415157 
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0.2413949 
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1.0000016114364 

nu 
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0.4500 

tauo 

= 

0.2415157 
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= 

0.2414191 

am 

= 

1.0000012754027 

nu 

= 

0.4500 

tauo 

= 

0.2415157 

tau 

= 

0.2414433 

am 

= 

1.0000009460213 

nu 

= 

0.4500 

tauo 

= 

0.2415157 

tau 

= 

0.2414674 

am 

= 

1.0000006236255 

nu 

= 

0.4500 

tauo 

= 

0.2415157 
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0.2414916 

am 
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1.0000003082242 
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0.4500 
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0.4500 

tauo 
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= 

0.2415399 
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nu 
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0.4500 

tauo 
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tau 

= 

0.2415640 

am 
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1.0000000000000 

nu 

= 

0.4500 

tauo 

= 

0.2415157 

tau 

= 

0.2415882 

am 

= 

1.0000000000000 

nu 
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0.4500 

tauo 
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0.2415157 
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0.2416123 
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1.0000000000000 

nu 
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0.4500 

tauo 
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0.3030615 

tau 
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0.5000 

tauo 

= 

0.3030615 
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nu 

= 

0.5500 

tauo 

= 

0.3732664 

tau 

= 

0.3730798 

am 

= 

1.0000000020424 

nu 

= 

0.5500 

tauo 

= 

0.3732664 

tau 

= 

0.3731171 

am 

= 

1.0000000010635 

nu 

= 

0.5500 

tauo 

= 

0.3732664 

tau 

= 

0.3731544 

am 

= 

1.0000000004564 

nu 

= 

0.5500 

tauo 

= 

0.3732664 

tau 

= 

0.3731918 

am 

= 

1.0000000001377 

nu 

= 

0.5500 

tauo 

= 

0.3732664 

tau 

= 

0.3732291 

am 

= 

1.0000000000175 

nu 

= 

0.5500 

tauo 

= 

0.3732664 

tau 

= 

0.3732664 

am 

= 

1.0000000000000 

nu 

= 

0.5500 

tauo 

= 

0.3732664 

tau 

= 

0.3733038 

am 

= 

1.0000000000000 

nu 

= 

0.5500 

tauo 

= 

0.3732664 

tau 

= 

0.3733411 

am 

= 

1.0000000000000 

nu 

= 

0.5500 

tauo 

= 

0.3732664 

tau 

= 

0.3733784 

am 

= 

1.0000000000000 

nu 

= 

0.5500 

tauo 

= 

0.3732664 

tau 

= 

0.3734157 

am 

= 

1.0000000000000 

nu 

= 

0.5500 

tauo 

= 

0.3732664 

tau 

= 

0.3734531 

am 

= 

1.0000000000000 

nu 

= 

0.6000 

tauo 

= 

0.4490661 

tau 

= 

0.4488415 

am 

= 

1.0000000004655 

nu 

= 

0.6000 

tauo 

= 

0.4490661 

tau 

= 

0.4488864 

am 

= 

1.0000000002389 

nu 

= 

0.6000 

tauo 

= 

0.4490661 

tau 

= 

0.4489313 

am 

= 

1.0000000001010 

nu 

= 

0.6000 

tauo 

= 

0.4490661 

tau 

= 

0.4489763 

am 

= 

1.0000000000300 

nu 

= 

0.6000 

tauo 

= 

0.4490661 

tau 

= 

0.4490212 

am 

= 

1.0000000000038 

nu 

= 

0.6000 

tauo 

= 

0.4490661 

tau 

= 

0.4490661 

am 

- 

1.0000000000000 

nu 

= 

0.6000 

tauo 

= 

0.4490661 

tau 

= 

0.4491110 

am 

= 

1.0000000000000 

nu 

= 

0.6000 

tauo 

= 

0.4490661 

tau 

= 

0.4491559 

am 

= 

1.0000000000000 

nu 

= 

0.6000 

tauo 

= 

0.4490661 

tau 

= 

0.4492008 

am 

= 

1.0000000000000 

nu 

= 

0.6000 

tauo 

= 

0.4490661 

tau 

= 

0.4492457 

am 

= 

1.0000000000000 

nu 

= 

0.6000 

tauo 

= 

0.4490661 

tau 

= 

0.4492906 

am 

= 

1.0000000000000 

nu 

= 

0.6500 

tauo 

= 

0.5277985 

tau 

= 

0.5275346 

am 

= 

1.0000000003183 

nu 

= 

0.6500 

tauo 

= 

0.5277985 

tau 

= 

0.5275874 

am 

= 

1.0000000001631 

nu 

= 

0.6500 

tauo 

= 

0.5277985 

tau 

= 

0.5276402 

am 

= 

1.0000000000689 

nu 

= 

0.6500 

tauo 

= 

0.5277985 

tau 

= 

0.5276930 

am 

= 

1.0000000000204 

nu 

= 

0.6500 

tauo 

= 

0.5277985 

tau 

= 

0.5277457 

am 

= 

1.0000000000026 

nu 

= 

0.6500 

tauo 

= 

0.5277985 

tau 

= 

0.5277985 

am 

= 

1.0000000000000 

nu 

= 

0.6500 

tauo 

= 

0.5277985 

tau 

= 

0.5278513 

am 

= 

1.0000000000000 

nu 

= 

0.6500 

tauo 

= 

0.5277985 

tau 

= 

0.5279041 

am 

= 

1.0000000000000 

nu 

= 

0.6500 

tauo 

= 

0.5277985 

tau 

= 

0.5279568 

am 

= 

1.0000000000000 

nu 

= 

0.6500 

tauo 

= 

0.5277985 

tau 

= 

0.5280096 

am 

= 

1.0000000000000 

nu 

= 

0.6500 

tauo 

= 

0.5277985 

tau 

= 

0.5280624 

am 

= 

1.0000000000000 

nu 

= 

0.7000 

tauo 

= 

0.6069004 

tau 

= 

0.6065970 

am 

= 

1.0000000003381 

nu 

= 

0.7000 

tauo 

= 

0.6069004 

tau 

= 

0.6066577 

am 

= 

1.0000000001732 

nu 

= 

0.7000 

tauo 

= 

0.6069004 

tau 

= 

0.6067184 

am 

= 

1.0000000000731 

nu 

= 

0.7000 

tauo 

= 

0.6069004 

tau 

= 

0.6067790 

am 

= 

1.0000000000216 

nu 

= 

0.7000 

tauo 

= 

0.6069004 

tau 

= 

0.6068397 

am 

= 

1.0000000000027 

nu 

= 

0.7000 

tauo 

= 

0.6069004 

tau 

= 

0.6069004 

am 

= 

1.0000000000000 

nu 

= 

0.7000 

tauo 

= 

0.6069004 

tau 

= 

0.6069611 

am 

= 

1.0000000000000 

nu 

= 

0.7000 

tauo 

= 

0.6069004 

tau 

= 

0.6070218 

am 

= 

1.0000000000000 

nu 

= 

0.7000 

tauo 

= 

0.6069004 

tau 

= 

0.6070825 

am 

= 

1.0000000000000 

nu 

= 

0.7000 

tauo 

= 

0.6069004 

tau 

= 

0.6071432 

am 

= 

1.0000000000000 

nu 

= 

0.7000 

tauo 

= 

0.6069004 

tau 

= 

0.6072039 

am 

= 

1.0000000000000 

nu 

= 

0.7500 

tauo 

= 

0.6840671 

tau 

= 

0.6837251 

am 

= 

1.0000000004786 

nu 

= 

0.7500 

tauo 

= 

0.6840671 

tau 

= 

0.6837935 

am 

= 

1.0000000002453 

nu 

= 

0.7500 

tauo 

= 

0.6840671 

tau 

= 

0.6838619 

am 

= 

1.0000000001036 

nu 

= 

0.7500 

tauo 

= 

0.6840671 

tau 

= 

0.6839303 

am 

= 

1.0000000000307 

nu 

= 

0.7500 

tauo 

= 

0.6840671 

tau 

= 

0.6839987 

am 

= 

1.0000000000038 

nu 

= 

0.7500 

tauo 

= 

0.6840671 

tau 

= 

0.6840671 

am 

= 

1.0000000000000 

nu 

= 

0.7500 

tauo 

= 

0.6840671 

tau 

= 

0.6841355 

am 

= 

1.0000000000000 

nu 

= 

0.7500 

tauo 

= 

0.6840671 

tau 

= 

0.6842039 

am 

= 

1.0000000000000 

nu 

= 

0.7500 

tauo 

= 

0.6840671 

tau 

= 

0.6842723 

am 

= 

1.0000000000000 

nu 

= 

0.7500 

tauo 

= 

0.6840671 

tau 

= 

0.6843407 

am 

= 

1.0000000000000 

nu 

= 

0.7500 

tauo 

= 

0.6840671 

tau 

= 

0.6844091 

am 

= 

1.0000000000000 

nu 

= 

0.8000 

tauo 

= 

0.7575481 

tau 

= 

0.7571694 

am 

= 

1.0000000008671 

nu 

= 

0.8000 

tauo 

= 

0.7575481 

tau 

= 

0.7572451 

am 

= 

1.0000000004446 

nu 

= 

0.8000 

tauo 

= 

0.7575481 

tau 

= 

0.7573209 

am 

= 

1.0000000001878 

nu 

= 

0.8000 

tauo 

= 

0.7575481 

tau 

= 

0.7573966 

am 

= 

1.0000000000557 

nu 

= 

0.8000 

tauo 

= 

0.7575481 

tau 

= 

0.7574724 

am 

= 

1.0000000000069 

nu 

= 

0.8000 

tauo 

= 

0.7575481 

tau 

= 

0.7575481 
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= 

1.0000000000000 
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= 

0.8000 

tauo 

= 

0.7575481 

tau 

= 

0.7576239 
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1.0000000000000 
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0.8000 

tauo 
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0.7575481 

tau 

= 

0.7576996 

am 

= 

1.0000000000000 
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nu 

= 

0.8000 

tauo 

= 

0.7575481 

tau 

= 

0.7577754 

am 

= 

1.0000000000000 

nu 

= 

0.8000 

tauo 

= 

0.7575481 

tau 

= 

0.7578512 

am 

= 

1.0000000000000 

nu 

= 

0.8000 

tauo 

= 

0.7575481 

tau 

= 

0.7579269 

am 

= 

1.0000000000000 

nu 

= 

0.8500 

tauo 

= 

0.8262315 

tau 

= 

0.8258184 

am 

= 

1.0000000020601 

nu 

= 

0.8500 

tauo 

= 

0.8262315 

tau 

= 

0.8259010 

am 

= 

1.0000000010570 

nu 

= 

0.8500 

tauo 

= 

0.8262315 

tau 

= 

0.8259836 

am 

= 

1.0000000004470 

nu 

= 

0.8500 

tauo 

= 

0.8262315 

tau 

= 

0.8260663 

am 

= 

1.0000000001328 

nu 

= 

0.8500 

tauo 

= 

0.8262315 

tau 

= 

0.8261489 

am 

= 

1.0000000000166 

nu 

= 

0.8500 

tauo 

= 

0.8262315 

tau 

= 

0.8262315 

am 

= 

1.0000000000000 

nu 

= 

0.8500 

tauo 

= 

0.8262315 

tau 

= 

0.8263141 

am 

= 

1.0000000000000 

nu 

= 

0.8500 

tauo 

= 

0.8262315 

tau 

= 

0.8263967 

am 

= 

1.0000000000000 

nu 

= 

0.8500 

tauo 

= 

0.8262315 

tau 

= 

0.8264794 

am 

= 

1.0000000000000 

nu 

= 

0.8500 

tauo 

= 

0.8262315 

tau 

= 

0.8265620 

am 

= 

1.0000000000000 

nu 

= 

0.8500 

tauo 

= 

0.8262315 

tau 

= 

0.8266446 

am 

= 

1.0000000000000 

nu 

= 

0.9000 

tauo 

= 

0.8895703 

tau 

= 

0.8891255 

am 

= 

1.0000000072633 

nu 

= 

0.9000 

tauo 

= 

0.8895703 

tau 

= 

0.8892145 

am 

= 

1.0000000037410 

nu 

= 

0.9000 

tauo 

= 

0.8895703 

tau 

= 

0.8893035 

am 

- 

1.0000000015878 

nu 

= 

0.9000 

tauo 

= 

0.8895703 

tau 
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